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3. Abstract •: ' '' ' OF £OOR QUALITY Sioux Falls, 50 b7 135 , 

■ This final report gives results on Landsat - 2 data received and evaluated 5 
scenes have been received in CCT Format and 65 scenes in film Format. These have been 
chosen to complete or extend time-series covering significant test areas. 

Data processing was mostly digital and different classification techniques were 
developped and comparatively evaluated with other available techniques. 

Output was in the form of colour or B/W maps using one or more scenes concur- 
rently. - ' • ■ 

Significant results were obtained in the following areas of investigation : - 

- coastal geomorphology, with possible applications in shellfish aquaculture or 
recreation devlopment ; 

- sedimentary transport, v>;ith identification of numerous sedimentary plumes 
CLoire estuary, Gironde estuary. Gulf of St-Malo, etc...] 

- methodological aspects of classification techniques. 
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1 - INTROnUCTIDN 


This report covers the whole of the FRALIT-B investigation period. 

A number of data processing techniques v/ere developed and applied 
to LAND3AT-2 CCT5. 

□n arrival, magnetic tapes were reformatted to our data hank standards 
(4 strips reassembled as a single scene file] and stored on direct access 
units. 

Classification of pixels in various spectral categories was followed by 
mapping and interpretation using ground truth. Spectal signatures were 
gathered on site using an EXOTECH radiometer and used in the training phase 
of the classification process. Emphasis was put on the mapping of coastal 
areas, extending seaward out to the limit of coastal sediment transports, 
and landward into salt marshes, dunes and polder environments. 

Comparisons were conducted in various test areas between LANDSAT data and 
other remotely sensed information (airplanes or other satellites) . Ground 
resolution as well as spectral band variations were assessed. 

Applications of remote sensing in coastal geomorphology were identified 
with particular reference to shell fish aquaculture and recreation ■ 
development . 



2 - TECHNIQUES 


2.1- COHPARATIVE ASSESSHENT OF DIGITAL TECHNIQUES 

The different classification methods used' in the FRALIT program were 
• applied to date, from a test site on the Fi'ench Atlantic coast. For 
critical evaluation of performance, this data was also processed by 
non-FRALIT methods CLARSY3, LARS JSC, lHAGE-100] . Before describing 
some of the results from this comparison, and the test site Itself, 
we will define some terms we shall use to characterize different 
techniques. 

2, 1.1 - Definitions 


Many different data processing techniques have been used in 
order to classify the spectral signature of remotely sensed 
data, into a given number of categories showing different 
terrain units or taxons [categories). 

In a classification process, operator participation fin the 
broadest sense) can occur at different phases and in various 
fields of competence. Therefore, we propose to use the follo- 
wing terminology : 

- . assistance : participation of the operator in data 
processing. This can take place at different phases. 

- pre-assistance : assistance occuring before the actual 
classificatinn process. 

- post -assistance ! assistance occuring after classification. 

- co-assistance : assistance carried out during the 
classification process. Pratice shows that the level 
of co-assistance is likely to increase with the speed 
and efficisncy of man-machine interaction. 

The experience required for operator. assistance can be drawn 
from different fields of competence. Therefore operator 
assistance can be termed spectral, taxonomic or geographical. 

- spectral assistance is based on the knowledge of the 
statistical characteristics of spectral signatures. 

The decision to subdivide a class having a bimodal 
spectral signature histogram is a case of spectral 
assistance, 

- taxonomic assistance requires prior knowledge of the 
physiographic units or taxons in a scene. An example 
is the decision to subdivide a forest class into 

.coniferous and deciduous. 
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- geographical assistance uses orior Knov^fledge of the area 
covered by a given type of landscape or of the limit 
between two farrain units. Geographical assistance, for 
example, is the decision to group into one class all the 
pixels belonging to a given ground 'area. 

Taxonomic and geograohical assistance are linked when they are 
based on the unambiguous association of a taxon and a clearly 
defined geographical area. 

Assistance is internal when it is based on the knowledge of the 
intrinsic oharacteristics of the data. The operator then partici- 
pates by assessing these characteristics in order to make dscisior 

Conversely, .assistance is external when data from other sources 
are used. This is the case when ground truth data is available 
and taken into account. 

The following table identifies the different processing techniques 
tried by FRALIT using the types of assistance as a classification 
criterion. 

2. 1.2 - Test zone description 

All the processing techniques used were tried out on the 
Fromentine test -zone for the following reasons [Fig 2.1. A] : 

1/ This area was imaged a number of times both by LANDSAT 1 and 
LANDSAT 2. Indeed, it had the lowest cloud cover of all the 
area studied in the FRALIT program. 

2/ It has a variety of coastal terrain units whose study is the 
subject the FRALIT program. Each of these units is large 
enough to be compatible with the spatial resolution of the 
remote sensors used. 

3/ It had been extensively studied by the FRALIT team. 

The Fromentine test site lies in the Vendee area and is defined 
by the following quadrangle : 

Let 46” 47’ N to 46° 57’ N 
and Long 2° 0 ’ W to 2° 20’ W 

It includes the following terrain-types : [Fig. 2.1.B] 

- Landward : Southern part of the Island of Noirmoutier 

: Breton Marsh 

• ^ Seaward ; Southern part of the bay of Bourgneuf 
: Atlantic Ocean [Bay of Biscay] 
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Fig. 2.1. A - Geographical location of the Fromontine test site 
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a - Land Envirnnmsnts 


a-1- Southern part of the Island of Moirmoutler 

This is basically a sand spit built against tertiary 
■ limestone out-crops. This spit includes coastal dunes with 
psammophyte cover locally mixed with maritime pine or 
Lambert’s cypress. 


East of dunes, a sandy marsh is divided up into small, 
variously cultivated stips. This marsh zone was extenreri 
by reclamation from the middle ages to the 19th century. ■ 
The latest dyKe lies north of the Gois ford : the Barbatrs 
polder holds nearly 300 acres and was,.hollowed out by the 
flood. The dyKe was patched round thV'feu'sti/ie -pf, this dip 
sc that a stagnant pond remains : thW ''’T''rdUJ de .Seh^a'^stopol" . 


.Sebastopol 


a-2- Breton Harsh 


The Breton Marsh is a marine alluvial accumulation 
[mostly clayl, protected by the dune covered bar of 
La Cote de Monts. The dunes are mostly built up along the 
sand bar to the West. Eastward, they lie across the bar 
and locally reach 23 metres height. They were planted 
with trees, mostly maritime pine, in the' nineteenth 
century. 


Housing is mostly developed where the dunes meet the marsh. 
Villages such as Fromentine and .INIotre Dame de Monts have 
fanned inland from the coast, with well-spaced' houses and 
numerous gardens . 

Land use in the marsh itself is very diversified : pastures 
in the east, with salt marshes [active or abandoned), sand 
or clay cultivated areas [Barre de Monts and La Crosniere) ' 
and occasional oyster beds [Bouin). Ancient polders have 
been invaded by the sea and partly reclaimed’ on the bay ■ 
of Bourgneuf coast. 

The are a is abundantly drained by narrow and winding 
canals, often linked to drainage ditches. 

< 

Winding roads also run through the niarsh. A four-lana 
highway crosses the south of the island of Moirmoutier, 
in line with the bridge linking the island to the continent 
near the narrowest point of the straits. 

b - Marine environment 


The marine parts of the test site are under the influence of a 
semi-diurnal tidevJhe tide drives a complex system of coastal 
currents which can reach 5 Knots in the Goulet de Fromentine. 
Tidal variations lead to the drying out of parts of the 
foreshore both in the southern bay of Bourgneuf and on the 
Atlantic coast. The Noirmoutier bridge marks the limit between 
these two environments. 



b“1 - Southern Bay of Bourgneuf 


From high to low tide marks, the foreshors can be 
subdivided into : 

. A polder, -Flooded by the sea in 1940 and never reclaim 
since : the Grande Prise de Fromentine. It has been 
largely silted up and the highest parts are covered by 
halophiles, mostly Suaeda Fruticosa. It is partly used 
as submersible oyster-beds. 

The Points aux Herbes de la Crosniere tidal marsh 
includes the following species from top to bottom : 
Suaeda Fruticosa on and near the dyke ; Pucclnella 
maritima and Sallcornia radicans, Obione portulacoides 
Suadea maritima and Aster Tripoliutn ; Salicornia 
herbacea and Spartina Stricta, 

, Flooded polders- with little or no halophile vegetation 
Beauvoir polder reclaimed in 1862 then enlarged in 16E 
and abandoned in the twentieth century. 

. Sandy mudflats behind the nineteenth century dykes. 
Geomorphologlcally active sand bars such as the Piau t 

. Mixed sand and mud areas on either side of the Gois- 
sedimentary ridge and its causeway. North of the. Gois 
these areas include oyster beds. 

. Tidal channels of varying depths above L-A.T. levels, 
running towards the straits in the southern part of 
the bay. To the North, the test-zone reaches the '-deep 
Fain French. 

b-2- Atlantic Ocean 

Qn the shores of the Atlantic Ocean, the following tide! 
areas can be distinguished : 

. Sand beaches on the Island of Noirmoutier, on the. Mon' 
coast and the shores of the strait. 

. Calcareo-siliceous limestone outcrops- : the la Fosse 
rocks, 

. Sand banks [submarine delta at Fromentine] for up to 
two miles seaward. These banks often include low 
frequency crosswise ridges with a changing morphology 

The water itself changes considerably with the tide and 
also with the weather. Three main types can be observec 
muddy waters in rough weather on a lowering tide in the 
Goulet de Fromentine i sand-laden waters- over the sub- 
marine delta on a rising tide s clear water off-shore o 
near the Cote de ilonts beaches. 

Sea state is variable : waves and swell are frequent in 
the Atlantic ocean, whereas the formation of foam takes 
place in the south part of the Bay of Bourgneuf VJhere 
v/aves are weak- 



“ 7 - 


2.2- CLASSIFICATION BY TAXONOmC PRE-AS5ISTAIMCE {FRACARTE } 
2. 2. 'I - Introduction 


The results presented here were obtained at the Ecole Nortfials 
Superieure (Montrouge) using the batch system FRACARTE tailored 
to the needs of the FRALIT team. 

Computation is done in batch processing on an IBM 370/1 B8. The 
test area is the Fromentine area as surveyed by LANDSAT2 on 
July 1975. 

. Computing time can be separated into two parts : 

- 12 seconds of CPU for legend construction during the 
pre-assistance phase which was conducted .by selection 
of training samples drawn from the entire population 
of 40 800 pixels. 

- 9.B seconds of CPU for the actual plotting file 
generation of a category map covering the 16 800 pixels 
test area. Off-line plotting was done in 45 minutes on 
a BENSON 122 ink-plotter. 

2. 2.2 - Software organization 

The system includes in two programs : ,FRALISIM and FRA.CARTE. 

FRALISIfl is used to identify the spectral signatures of the 
Know taxons located in the area of interest. This category 
determination is drav^'n from the selection of sufficiently 
homogeneous training areas in the image to be processed. 

When the geographical area to be catagorized is too small, 
it is necessary to select training samples outside it in 
order to draw conclusions from homogeneous taxons covering 
a statistically significant number of pixels. The spectral 
signature determination is done on gray scale prints or plain 
listing outputs. These listings help in the determination in 
each spectral band of the radiance intervals 'which bound the 
multispectral signatures of each taxon. 

Once this determination is complete the spectral characteristics 
of each category are input on punched cards to the FRACARTE 
program. This module generates a magnetic tape plotting file 
which Is used to drive the off-line BENSON plotter. Any orojectior 
and scales can be selected. In the example presented the operator 
selected a Lambert2 projection and 1:50 DDO scale. 

2. 2.3 - Theoretical backgroung 


All the philosophy of this system relies on a high level of 
pre-assistance used for the determination of radiance intervals 
bounding the spectral characteristics of each recognized 
category CtaxonO . 



This detarmination can be guided using two types of approaches : 

a - Assuming the taxons belong to spectrally well defined and 
homogeneous categories. Their geographical boundaries as 
well as their spectral limits are clearly defined because 
the spectral gradients at their boundaries are strong. In 
LANDSAT data this type of signature is more frequent on 
land than at sea and is very often associated with man-made 
environments rather than natural ones. However it is possibl 
to find moderate gradients at sea [Fig.- 2. 2. A) especially 
in estuarine or coastal areas. 

In this case, the taxonomic pre-assistance will produce 
unambiguous spectral bounds for each category. It is seldom 
necessary with LANDSAT to use the entire set of spectral 
bands, 2 spectral bands are sufficient in most cases. 

b - The taxons cover continuously varying environments which 
are really obtained by mixing in various proportions a 
limited number of categories. The limits are in this case • 
rather fuzzy and variations of signatures do not generate 
significant gradients. 

We propose to call them spectral continuums. They can be fot 

- In coastal waters, especially in the Inside of bays 

- In the intertidal environments 

In this case, the taxonomic pre-assistance can usefully be 
followed by a spectral pre-assistance working on the statis- 
tics of data. For example, types of waters defined by 
spectral intervals can be mapped by increasing turbidity as 
this parameter is correlated with spectral bands HSS4 and 
HSS5. 

The weak points of this approach are : 

- its dependance on careful choice of training areas 
in which categories are Known for their homogeneity 
by collecting appropriate ground truth. ■ 

~ its reliance on the "a priori" knowledge and skills 
of the operator. 

But. it is also very useful because of : 

- its simplicity of use 

- its small computing times 

- and its results in the case of spectral continuums j 
classification because the produced maps are equivaler 
to photographic isodensity slicing. 



Fig. 2. 2, A - Example of coastal water colot’ gradient 
At ebb tide, the plume of turbid waters 
coming out the Maumusson strait expands 
over clear ocean waters [£.2331 - 10110, 

19 Dec 75, line 156K3, 

Values in MSS7 werei'^not drawn because always 
equal to zero or one. 
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2. 2.4 - An example of processing 

The first step of the processing of the Fromentine ’area leads tc 
the determination of class spectral intervals drawn from 
■homogeneous training samples. The listing output of FRASILIM 
gives a list of values for the four HSS spectral 
bands. Each individual symbol is assigned to a specific radlancs 
value (Fig. 2.2.B1 . 

From this, listing the spectral signatures of 6 categories were 
identified : 



- water 

- tidal flat's 

- sand and built-up areas 

*- dense forest 

r cultures, meadows, vegetated sand dunes and open forest 

- dispersed cultures and settlements with gardens 

• These definitions are assigned to each signature by taxonomic 
assistance'. Two categories from this set can be interpreted as 
spectral continuums : the water and tidal flat environments. 
■ThBoretical considerations' on spectral signatures or ground- 
truth radiometric measurements done with an EXOTECH radiometer 
lead to the subdivision of the first environment into 5 sub- 
categories of water ranked by their radiance values in spectral, 
band MSS5 : 


- MSS 5 


. 23 

class ;1 

- MSS 5 


24-27 

class 2 

- MSS 5 


28-31 

- class 3 

- MSS5 ' 


32-35 

class 4 

- MSS5 


3B 

class 5 


'These waters are ranked by increasing turbidity level. The tide 
flat category extends its signature in (1537 from radiance value 
B to 12. It is divided into two sub-classes using the intervals 
[6,9j. and [^10,12j . This gives us a separation between humid 
slikkes and tidal flats with lower water content. 

The previous taxonomic pre-assistance- has been refined by 
spectral interpretation of the spectral continue . The subsequen' 
intervals are input to FRACARTE which outputs a 1:50 000 scale ■ 
map. (Fig. 2'.2.B]. The legend of this document can be described 
as follows. 

1 = Clear marine waters. They extend off the Monts coast 

(lower left part of map) . 

2 = Slighly turbid marine waters in Bourgneuf bay and Nort! 

of the Gois ford Cupper left part of map) as well as Oi 
the outer edge of the Fromentir|e submarine delta. 
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3 = Averags turbid marine waters. 

4 = Highly turbid marine waters. 

5 = Marine waters with very high sediment content in 

suspension often mixea with slightly covered tidal flat 

G = Very humid slikkes. 

7 = Tidal flats with low water content. 

8 = Maritime pine forest and salt marshes (Monts forest 

and salt marshes on the rim of the Bourgnsuf ■ bay noest] 

0 = Cultures, open forest, meadow and vegetated sand dunes. 

10 = Dispersed cultures and garden settlements. 

11 = Sand and built-up areas. 

Objects of different natures such ss the maritime pine forest of 
Monts and the salt marshes .of the abandoned Grande Prise golders 
may have the same spectral signature. In fact, three .pixels 
belonging to these areas display identical quadruplet signatures 
23, 24, 37 and 16 respectively in spectral bands MSS4, 5, 6 and 
(cf. quadruplets NM=G in Fig. 2.2. C). Resolving such ambiguities 
may involve the use of time varying signatures (diachronic) or 
spatial signatures [texture]. 

Results can be further refined on land but mlsclasslfication 
may increase due the complexity of the agricultural pattern of 
the Marais de Monts and Noirmoutier island. We decided to i^eep 
this level of category division because it is well correlated 
with physiographic units and broad categories of land use. 
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Fig. 2.2.B - An example of classification by taxonomic pre-assistance 
Map of the Fromentine area at 1/50 000 scale generated 
from LANDSAT2 data (E2183-1Q161 - 29 July 1975] 
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Fig. 2.2.C “ Part of a FRALISIH output listing for the 

Fromentine area C29 Juillet 197B] . Each quadruplet 
gives the value taken by one pixel in the four 
spectral bands MSS4 to HSS7 . 

The coding of individual values is hexadecimal. 

The determination of category CtaxonD signature 
is taken from such a document during taxonomic . 
pre-assistance. The quadruplets Nf1=G (values 23, 24, 
37, IB] have been circled out. They are associated 
to the Monts forest [last line of display] or the 
Grande prise salt marshes [first line of display] 
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2.3- CLASSIFICATION BY SPECTRAL P RE -ASSISTANCE CAGREG ] 
2-‘.'3.1 - Introduction 


Data processing resulrs djscuBsed here, were obtained using 
the AGREG program which is part of the CLAMS classification 
system. This system v;as developped by Ecole des Mines under 
contract with the Centre National d ’Etudes Spatiales CCNES] . 
It is implemented in batch processing at the CNES computing 
center in Toulouse, using a CDC 7600. Overall CPU computing 
time for treatement of the Fromentine area was about 1 minute 
- and 20 seconds. 

Graphic output related to the classification results is 
generated as printed lists and graphs drawion black and White 
micro.fi Im using a CALCOflP 1670 COM system. 


2. 3.2 - Softvvjare organization 




The use of the AGREG algorithm causes the execution of several 
modules integrated in the CLAMS system. They perform the 
following actions : ■ 

- supervised clustering per se CAGREG] 

- image display of classification results 

- spectral, statistics 

"A priori" knowledge by the operator of the number of classes 
and of their spectral statistics [average and RMS) is needed in 
order to run AGREG. These' characteristics can be acquired using 
various means of information ; 

- a spectral data bank built by collecting typical spectra: 
signatures. In this case the statistics are estimated by 
the variability expected from natural environments, 

- a previous classification applied to a sample of the 
original image using a limited "number of pixels. The 
reduction in the number of objects to be classified 
would in this case be linked to the use of more 
sophisticated methods than AGREG thus more costly in 
computing time. 

- a previous classification carried out on a different are; 
or even on data acquired for the area of interest but at 
a different date. The philosophy is in this last case to 
extend 'in time end in space the results obtained on an 
intensive test site where ground thruth was thoroughly 
collected. 


Once clustering has been performed around cluster centers 
initialized by spectral pre-assistance, category display and 
spectral statistics help in assessing the quality of the 
subsequent mapping. 
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2 . 3.3 - Theoretical background 


The supervised clustering performed by AGREG belongs to the 
general .family of algorithms in which spectral radiances are 
assigned to classes v;hich are ”a priori” defined by their 
statistical characteristics. 

Several practical approaches are found in order to solve this 
type of problem. They are generally of three types : 

- assignment by computation of the conditional probability 
that a spectral signature belongs to a specific class 

- assignment by reference to a decision table 

“ assignment by computation of the distance separating a 
spectral signature and a specific class. 

a - Assignment by conditional probability 


The first type of approach relies cn the theory of 
probability. The basic information used in this case is the 
joint radiance histogram computed for the entire set of 
spectral bands..' 

This basic set of data can be information compressed at variou: 
degrees depending on the available comiputing time and core 
metnory space. This compression usually involves the use of 
probability distribution modelling of the spectral radiance 
values repartition of pixels known to belong to specific 
. classes. This training is done by ”a priori” gathering of 
ground truth. 

The model may use a limited number of parameters. This is 
for example the case of Gaussian hypothesis used in the 
LARSYS system. Illustration of this type of parametric 
approach can be found in Fig 2. 3. A in the case where two 
spectral bands are used. A more precise but more complex 
solution can be found by the use of general types of 
distributions as in the general class of non parametric 
methods . ' . . 

b - Assignment using a decision table 


The basic information used Is the same as in the preceeding 
methodj namely the Joint radiance ; histogram. 

But in this case no effort Is done to model this distribution.- 
Data compression is only performed in the storing phase in 
order to built a multidimensional decision table [Fig 2. 3. A]. 

The assignment' is in this case very quick but core memory 
used is often in the order of several tens of thousands of 
bytes. Memory reduction is in this case enterily depencent 
on computing skills developed by trained personel. 



ZSSIAI 





MSS4 


Fig. 2. 3. A - Two types of assignment methods using 

the joint radiance histogram of spectral 
signatures 
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C - Assignement using a distance criterion 

The degree of similarity between a spectrum and a class is 
computed in this case in place of the conditlcnal probabi- 
lity of a spectrum to belong to a class. The solution used 
.in general is to compute the Hahalanobis distance which csss 
for each class C., its average signature m. and the associa- 
ted variance “Covariance matrix ^ 

D„Cz,C J = fz-mJ 

n i XI 1 

where : E, = inverse of matrix E. 

X X 

This formula defines a quadratic form associated to the 
nahalanobis metric. 


The AGREG program designers prefer a distance which uses 
the average signature m_, and the vector a.' which set of 
components a are the spectral variances^class C. for the 
spectral band k. It is defined as follows : ^ 


□ Cz,Cj 

X 


n 

= E 
K=1 




/ a 


iK 


The use of this distance gives classification results which 
are equivalent in quality to the use of D|^. However it is 
more simple to compute than the last distance, leading to 
a saving in computer time. 


A reject class (unclassified points] is built using a crite- 
rion equivalent to the mean divergence. 

nethods of assignement using similarity indexes do save the 
computing time used in fitting specific models to the 
radiance joint distribution histogram but this saving is 
often lost when the assignment phase involves a large 
number of pixels. 


2. 3.4 - Data processing example 

a - ’’Apriorl” class statistics 

Decision was made to categorize the Fromentine area into 
six "a priori” classes [Fig 2.3.B] whose statistical charac- 
•terlstics ware drawn from a detailed classification computed 
on a regular sampling of the studied area. 

The six classes are related to land use themes and do not 
cover sea water or tidal flats. 

b - Classification performance evaluation 


A solution to the "a posteriori" evaluation of the 
supervised clustering performance may be found in the 
comparaison of the position of means and variances of 
"a priori” classes in relation to the "a posteriori” 
spectral signature histograms computed in each specxral band. 



MSS7 



Fig. 2.3.8 


1 2 


. B 

spectral characteristics of the classes 
involved In the supervised clustering phase 

A - Spectral means and variances for each 
group for spectral bands iiSS4 and MSS7 

B - Average spectral signatures 

>.^.>...3 — _C 
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This test was performed on p3.g 2.3.C comparing the MSS4 
histograms of groups 1 and 2. The hatched part In sacn 
histogram is related to the radiance interval [m-a, m+a] 
computed from the "a priori" spectral characteristics 
of the associated class. 

A good fit is noticeable for group 2 between posteriori" 
histogram and ‘"a priori" spectral characteristics. 

This is not tne case for group 1 for which the histogram he 
a trend toward blmodallty . This property was not taken into 
account in the similarity distance formula used in AGREG 
because it is only computed from spectral mean and variance 
A probability assignment method, using a non-parametric 
approach, might have given a better result to the expense 
of more computing time. 

c - Results display 


The cartography presented on Fig 2.3. D is an enlargement 
of a display ^generated on microfilm. 

The reject class is represented in white. It contains the 
pixels which spectral radiances were . found to be too dlffs 
rant from those of the "a priori" classes'. 

The taxonomic significance of the classes was found to be 
in good relationship with existing land use patterns. 
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Fig. 2.3.C - Comparison of ”a postsriori” spectral 
characteristics of groups 1 and 2 with 
the "a priori" parameters used in the 
clustering phase 
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Enlargenfient of AGREG classification display 
generated on a COIl mlcrcfilm device 
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2.4- CLASSIFICATION BY SPECTRAL CDASSISTAMCE AND TAXONOMIC 
POST ASSISTANCE (FRACAMT 


2. 4.1 - Introduction 

The data analysis results presented here are taken from the 
Fromentine area, overpassed by LANDSAT1 on July , 1973 
(Ref. number: 1353-10244]. The equiment used is a CII IRIS-80 
computer linked in the time-sharing mode to a Tektronix graphic 
display terminal. The graphs and maps produced on this interac- 
tive console can be hard-copied on paper using ‘a Versatec 
printer-plotter. 

The overall computing time for a classification involving 
40 000 pixels is about 20 minutes. 


2. 4.2 - Software organization 

'The FRACAM system has been developped at the Ecole des Nines 
(Paris). It is organized in four separate modules, whose 
functions are : . 






- data display and spectral statistics (-FRAC1), 

- correspondence- analysis '(FRAC2), 

- ascending' hierarchical classification (FRAC3), 

- adaptatiVB clustering (FRAC4) . 

The FRAC1 program has various display capabilities whereby 
images can be presented as : 

• s 

- black and white binary maps related to a specific interve 
of radiance values in a- given spectral band or to a 
spectral class produced fay a previous analysis. 

- gray-scale images, with a limited number. of gray scales 
(4 levels at maximum), fay modulating the intensity of the 

. electron beam impinging on the storage tube. 

- as an infographic map using the same logic as the FRACARTj 
program. 


This program produces graphs related to the following statis- 
tical results : 


- radiance values histogram, in a specific spectral band, 
computed on the entire image or restricted to a spectral 
class produced by a previous analysis. 

- average and RNS .signatures of pixels associated to a 
spectral class or to an interval of radiance values in 
a given spectral band. 

- bidimenslonal histograms obtained by the comparison of 
two spectral hands (f1SS4 and NSS7, for example) . These ce 
be computed for a complete image or can be restricted to 
spectral class. 
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- a display of . spectral class signaturas in the plane of 
two spectral bands. 

The maps produced by FRAC1 are geometrically corrected using 
LAND3AT orbit parameters and scanning geometry specifications. 

The hard-copies generated on the printer-plotter match the 
Lanibert 2 projection system at any scale between 1t1000 000 
end 1 :10 000. , - 

Another function of the FRAC1 program is to- 'prepare the .coding 
of data, either by slicing the histograms In- isopooulatlon 
intervals, or by adapting these intervals to the spectral 
contrasts the operator is able to assess in the image before 
any classification (opposition between water and soil, for 
example] . The first step in any pattern recognition problem 
is to select an appropriate feature space by data coding end 
also. In most cases data compression. This is done in the FRAC2 
program. Feature space selection uses a two step procedure where 
data ceding is followed by data compression. A binary coding of 
the radiance values using the previously determined histogram 
intervals is followed by correspondence analysis. This multidi- 
mensional statistical method similar to factor analysis enabiss 
the operator to select a sub-space of limited dimensions (5 typi- 
cally] in which hierarchical classification and clustering will 
be performed. This analysis can be carried out for. ^the-.entire 
image or can be restricted to a spectral class produ'cedi-by a 
previous analysis. 

The FRAC3 program helps in initiating the clustering of radiance 
spectra in groups of similar shapes. The number of groups that 
one would expect to find in the data set is objectively determined 
as 'a result of an ascending hierarchical classification of the 
spectral intervals used in the data coding performed in FRAC2. 

Dialogue between the computer and the operator is required in 'this 
phase, because the.determination of the expected number of groups 
involves the Interactive "pruning” of a dendogram resulting from 
the hierarchical classification. Along with the -number of groups 
the position of their centers of mass in the .feature space is 
estimated as well. 

FRAC3 output is used in order to Initialize the adaptative 
clustering performed by FRAC4 . It uses an iterative procedure 
of the same principle as the ISODATA clustering algorithm. 

Within each iteration every radiance spectrum is assigned to 
a specific group whose center of mass is updated at the end of 
the assignment phase. Various convergence cri-teria are used 
such as an upper bound on overall center displacements or total 
number of class assignments changes between e.ach iteration. 
Operator action is limited to the choice of the maximum number 
of iterationB in the case vjhere convergence does not occur 
using built in criteria. 

The FRACAfl system is organized so as to enable the operator to 
perform a layered classification by sequentially refining rhs 
spectral groups [spec.tral co-assistancel in relation to their 
taxonomic significance- -(^taxonomic post-assistance] . 
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2. 4^3 - Theoretical background 
a - Methodology 

The relationship existing between an object located on the 
surface of the earth and the shape of its radiance spectrum as 
tnesured by LANDSAT is the basic principle of remote sensing. 

It lends itself to automatic cartography through the use of 
numerical procedures which by mathematical nature belong to' - 
pattern recognition techniques. 

The recognition of spectral shapes is solved by two types of 
methods. Segmentation of the set of spectral radiances is 
based, in the supervised approach, on the "a priori" knowledge 
of the nature and location of typical objects existing in a 
given test area. In the unsupervised approach, clustering using 
objective criteria. 



The classification method used in FRACAM belongs to the latter 
type. It is built around an adaptative clustering algorithm 
similar to the ISDDATA method. 

As the number of objeots to be recognized can be rather large 
(8 X 10^ for a complete LANDSAT image), special care has been 
taken to reduce computing time. For an ISODATA type algorithm 
this mainly requires the minimization of the number of iteratic 
necessary to reach convergence. This problem is equivalent in 
some respects to the search for Initial centers of mass as clos 
as possible to their final positions after convergence. 

Correspondence analysis of J.P. BENZECRI as a multidimensional 
method for data structure analysis provides an answer to the 
search for reliable centers of mass. It results in the selectio 
of a sub-space of limited dimensionality (feature space] in 
which it is possible to undertake a detailed study of the 
similarities existing between the spectral bands or their coded 
equivalents. The groupings thus determined are used in the 
.initialization phase of the clustering algorithm. Initial 
centers of mass are computed using the simultaneous representati 
formula of correspondence analysis which builds a straightforwa 
link between the factor space of the spectral bands and th'e Tac 
space of the spectral signatures. 

Using this general framework, the FRACAM system is logically 
organized in three steps ; 

- use of correspondence analysis in an objective study of 
the data set structure followed by the selection cf a 
factor sub-space which is retained as feature space. 

This data compression very often filters instrument noise 

- determination on the filtered data of an ascending hierar 
chical classification of the spectral bands or their code 
equivalents. This similarity study leads to the recogniti 
of associations between spectral bands and spectral shape 

- construction of spectrally homogeneous groups by adaptatl 
clustering where the initial positions of the centers of 
mass in the feature space are calculated from the results 
of the previous phase. 
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b - D ata coding 

In order to have more flexibility in the determination of 
spectral band groups, it is preferable to transform the 
originally measured spectral radiances. This transformation 
is called complete disjunctive coding or isopopulation coding. 

Such a code is generated by the following steps : computation 
of the radiance histogram of 'each spectral band j for the 
area of study determination of n, isopo'oulatlon radiance 
intervals for each spectral band 3 and replacement in the 
original data set of each measured radiance spectrum by its 
complete disjunctive code. 

Each n dimensional spectrum Cn = 4 for LANDSAT) is thus ' 
replaced by a K-dlmensional vector, with : 

K. =^|ri^ i j = 1 

Each component of this vector associated to a given radiance 
interval takes a binary value 1 or 0 depending wether the 
corresponding spectrum radiance value in the spectral band j 
belongs to the given radiance interval or not. 

The principle of this coding Is illustrated in Fig 2. 4. A 
•v.ihen the four spectral bands of LANQSAT are snlit into four 
isopopulation intervals. 

The use of this procedure was justified by en increase in 
overall classification performance, specifically regarding 
the number and spectral definition-of classes. 

c - Correspondence analysis 



After transforming the data set. it is interesting to study 
its structure with particular reference to the similarity 
existing between the K coded equivalents to the n initial 
spectral bands. 

We will call transformed band c. , the binary instrument 
whose output is 1 or 0, wether'^^he observed object radiance 
in band j falls or not in the Isopopulation radiance interval - 
of rank a. 

In the example of Fig. 2.4.B the transformed band a 

very simple incremental instrument whith an output of 1 v^hen 
the measured radiance in band 4 fMSS7) falls between intensi- 
ties 4 and 1B or □ for any radiance value outside this intervs 

The total number of transformed bands is K. (K='1B in Fig. 2.4.,-’ 
they will be indexed by k in our formulas as long as the 
explicit identification of the relationship Detween spectral 
bands and transformed bands is not required. 
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The study of the similarity betv\'een two transformed bands K 
and k' in done by using the following algebiaic structure. 
The transformed bands are considered as points of an I 
dimensional vector space in which I is the number of 
pixels to be classified. The origin of this vector space 
is the average radiance computed in each spectral band 
for the entire set of pixels. 

Fig. 2.4.B gives an example of the data table associated 
to the situation shown in Fig. 2. 4. A. If one seeks a 
mechanical equivalent to this data setj it can be th-ought 
of as a ’’cloud” of k transformed bands scattered in the 
vector space E^.. 


Each point k belonging to this ’’cloud” in is given a 
weight p which is the global energy received by the 
incremental band K. The metric of E^ is defined by a 
similarity distance called distance (Chi square distance 
defined as follows i 


with. : 



Vm Z 




1 





pik'‘p 


k 


- n = number of actual spectral bands (n=4 for 
LANDS AT] 


This distance is equal to zero vi/hen the binary images 
associated with k and k' are identical. Theoratically the 
.computation of the distance existing between each couple 
(k,k’] is sufficient for the study of their similarity. 

2 

Unfortunately the computation of in Ej is long and 

tedious because of the very large number of dimensions 
of this vector space. Typically several billions for a 
LANDSAT image. 


Fortunately, correspondence analysis enables us to select 
a feature space F of limited dimensionality (ra=5 for examp 
F is a sub-space of which is choosen so that in F the 
dxstance between the projections of the transformed bands k 
and k’ is a good approximation of the true distance d , , 
between the two points E^. The base vectors of F are^the 
eigenvectors of the "cloud" built v;ith the K transformed 
bands weighted by the . The associated eigenvalues A 
are ranked by decreasing^order : 


X 


1 


> X. 


>. . . 


> X 

tn 


Each of them can be interpreted as representing the percen- 
tage of information Cor inertia) born by the associated 
eigenvector . These eigenvectors are the coordinate axes 
of F which is often called factor space. It is in F that 
the other steps of FRACAM take place. ^ 
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d - Typology of transforaed bands 


In order to find reliable initial conditions for the 
adaptativ/e clustering of spectra in F j it is necessary to 
have an estimation both of the number of groups they form 
and of the positions of the corresponding -centers of mass. 

Our philosophy is to consider that groups of similar 
transformed bands are most likely' to characterize specific 
types of spectral shapes (or signatures] .• This is why 
groupings of transformed bands are sought for in F^. 

Using the simultaneous representation formula- of correspon- 
dence analysis we will then be able to associate to the 
center of mass of each group of transformed bands a point 
in the spectral space which is a good approximation of the 
"true" center of mass of the corresponding group of similar 
spectra. 

Before performing this spectral band to spectral signature 

transition, a detailed study of the geometric structure of 

the "cloud" built by the projections of the transformed 

bands is undertaken. This is done by using an ascending 

hierarchical classification algorithm in F . 

m 

The construction of a hierarchy on the set on K projections 
is dona by gradually building a series of mutually inclusive 
classes fitted one into another. A dendogram representation 
Is given to this series (Fig. 2.4. C] in which each node is 
indexed by a number proportional to the diameter in F of 
the corresponding class. The ascending character of tins 
selected algorithm comes from upward construction of the 
dendogram (tree) , the classes being formed by sequencially 
grouping points or groups of points which are close to one 
another. 

The determination of the characteristic groupings between 
similar bands is done by "cutting” the dendogram at a given 
diameter by the operator thus defining an upper bound to 
the diameter of the classes being constructed. 

In the case of Fig. 2.4..C three groups of similar transformed 
bands are selected as' : 


- group 

1 

■ '^44’ 

S2' 

S4- 

S4' 

- group 

2 


Ss' 



— group 

3 

• S3' 

S2’ 

S2' 

Si ' 



S2' 

S4’ 

S3' 

S3- 


A center of mass (center of gravity) is associated to each 
of those groups as a points in F . The spatial distribution 
of the groups of transformed ban9 "points" corresponding to 
the previous determination is presented in Fig. 2.4.D. We 
have limited ourvelves to a plans representation restricted 
to the first and second factors (m=2) . 
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Fig. 2.4.D - Display of a hierarchy of coded spectra-l 
bands generated from the results of 
Fig. 2.4.C. 

The representation is done in the plane 
of factors 1 and 2'. 
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e - Clustering of spectra 


One of the more important properties of correspondence 
analysis is to solve at the same time the twin problem 
of structure analysis of similarities both between spectral 
bands and between spectral signatures- In fact, the set of 
eigenvalues | a = 1,mj- which is used in order to built 
F can also be%sed to built a feature space of the 
measured spectra with tne same number of dimension as F . 
This sub-space G is constructed from F by simple multipli- 
cation of the base vectors components by the square roots 

of the associated eigenvalues X «• 

GC 

Clustering taKes place in G using an adaptative centers 
of mass algorithm which borrows some of its properties from 
the K.-mBans and the ISDDATA algorithms. Initialization takes 
■place by specification of the number of groups and initial 
positions of the cluster centers. 

In the present case the number of spectral groups is equal 
to the number of groups of similar bands as determined by 
the previous phase. The initial cluster centers are the 
projections in G of the centers of mass of the groups of 
similar bands determined ' in F . Projection is easily done 
by multiplying each of the m coRlponents of the centers of 
mass by the corresponding eigenvalue 

The clustering algorithm is stoped, and radiance spectra 
are assigned to a class by using a minimal distance criterioj 
in G . Namely, each spectra is assigned to the closest cluste| 
center. After each iteration, the cluster centers are update' 
with regard to the effective population of each class. 

The philosophy of this approach is 'that after a certain 
number of iterations, cluster centers and cluster cardinals 
are sufficiently stable to give a reliable classification. 
Various stability criteria are computed at the end of each | 
iteration in order to assess the quality of the classif icatij 
The criteria are : the ratio between intra-class variance an 
inter-class ■ variance, the overall displacement of cluster 
centers between their previous and updated positions and the 
number of class assignment changes from one iteration to 
the next. 

At the end of previous algorithm, various results are avails 
First, a cartography in which each pixel is given a symbol I 
uniquely associated to the class to which the associated 
radiance spectrum belongs. Statistical results are also 
produced for each class in the form of average and R.M.S. 
radiance values for eaoh spectral band. 

Experience has shown that the initial cluster centers, as 
determinated from the transformed bands typology are very 
reliable. This has enabled us to reduce the number of itefo' 
tions to 2 our 3 thus obtaining a considerable saving in 
overall computing time. 
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f - An adaptatlve descending classification 

In the pratical use of unsupervised clustering algorithms, 
two approaches are possible. 

In the first case, a- large number of clusters is determined, 
-With the risk, of having to followup the- clustering phase 
by a similarity study of the groups in order to reduce their 
number. In the second case, the clustering algorithm procseds 
by successive dichotomy until the desired number of classes 
is reached. 

The first type of heuristics is used in the LARSYS system, 
in which the first step leads to the determination of a 
•large number of- potential classes taken from the training 
area. The study, of the similarities existing between those 
■ groups of limited cardinals, enables to regroup them follo- 
wing a typical phase of taxonomic pre-assistance. This method 
uses an ascending approach whereby the complexity of a weak 
data compression Clarge number of groups) is reduced by the 
relative simplicity of a taxonomy often guided by typical 
application themes (oceanography, agriculture, geology ...)'., 

The algorithm used in ISOCLS belongs to the second type of 
approach, it proceeds by successive dichotomy between' each 
iteration. This is a descending approach going from the 
general to the particular. 

The structure of the FRACAM system looks well suited to the 
second type of heuristics. The ''cutting" phase of the 
hierarchical dendogram enables the operator to control the 
number of groups while performing a polytomy which is more 
complex than ISOCLS. This descending classification process 
is also made adaptative because at each downward move, 
correspondence analysis and hierarchical classification are 
computed again in order to adapt their results to the statis- 
tical properties of the group which the operator wants to 
refine. 

Several criteria for assessing the need to pursue the subdivi- 
sion of a group can be used. We propose to rely on the shape 
of the radiance histogram in MSS7. The subdivision is pursued 
if the group histogram is multimodal or does not fit a simple 
probability distributicn model (Gauss, unifoim or lognormal 
distributions): For example, we have reproduced in Fig. 2.4.E 
the spectral distributions of the class of pixels obtained 
after a first classification of the Fromentlne area on July 11, 
1973. Six groups were determined, the use of above criterion 

leads us to the following conclusions : 

' 

- groups 4, 5 and 6 do not need to be subdivided 
-•groups 1, 2 and 3 are on the contrary eligible for 
futher processing 

Going from' divisions to subdivisions, the FRACAM system Is 
able to perform an adaptative descending classification guided, 
by a spectral coasslstance taylored to the group spectral 
statistics. 



Fig. 2.4.E “ MS37 spectral band histograms for the 

six groups obtained after the first phase 
of treatment precseding spectral co-assistance 
The unnumbered plain line is the nsS7 histogram 
for the entire image. 
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2. 4.4 - An example of processing 


The processing of the Frcmentine area surveyed on July llj “\B73 
was done in two steps. In the first step, complete disj.unctiVB 
coding of nsS4 to WSS6 in eight intervals and MSS7 in seven 
intervals led to the determination of six groups. Spectral 
coassistancs was performed on the band fiSS 7- histograms (Fig. 2. 4. E] 
thus resulting in the further subdivision of groups 1, 2 and 3. 

The adaptative descending classification was pursed for these 
groups, by first generating a new complete disjunctive code 
with four intervals in each spectral band. This code was adaptacec 
to each group's statistics seperately. Finally, eleven groups 
were determined by the classification and submitted to interpre- 
tation in the taxonomic post-assistance phase. 

a - The results of the first phase 

Regarding correspondence analysis we present on Fig. 2.4. F 
the projections cf the transformed bands in the plane of the 
first and second factors. In this representation, each spec- 
tral band (MSS4 to MSS7) was assigned a letter (A to D) and 
a number corresponding to a complete disjunctive interval. 

For example B7 is the the transformed band related. to the 
seventh^ spectral interval determinated from the MSS5- histogram 

Each factor can he interpreted as follows : 

- on the first axis, the general structure corresponds 

to a transition from low radiance values to high radianc 
values in spectral bands MSS5, flSSB and MSS7. 

- on the second axis, spectral .band HSS4 is organised 
from low to high values except for the near center 
position of A1 (radiances from 2B to 32 in f1SS4] . 

The two dimensional plane representation enables us to 
refine this interpretation. A strong correlation is observed 
between NSS6 and MSS7, The parallel structures drawn by these 
spectral bands are organised in two arcs of opposite curvature 

- a first arc spreads over factor 1 going from C , D to 
^ 4 ' °3 

- a second arc spreads over factor 2 going from C^, 
to Cq, 

This general organization is due to a broad dichotomy in the 
data set between two types of spectral signatures correspon- 
ding, in the surveyed area, to the opposition between water 
themes (sea, channels, tidal flats) and land themes (cultui'es, 
forest, marshes] . 

The ascending hierarchical classification has led to the 
determination of six groups of transformed bands v;hich are 
oultined in Fig. 2.4.F. The resulting, spectral clusters were 
displayed one groups at a time as binary maps. One can find 
an example of a such a map in Fig. 2.4.G. 
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Fig. 2.4.F - Display' of first phase results in factor 
plane 1-2. 

Groups^- of coded equivalent spectral bands 
were determined by hierarchical ascending 
classification . 



3B 


0 ^ 




This display of group B was obtained as a Versatec printer- 
plotter hardcopy of a Tektronix graphic display. 

The spectral coasslstance phase involves the active particl 
pation of the operator in assessing the spectral individual 
□f each group. This is done using spectral signature statis 
for each group : radiance values histogram, average spectra 
signature or even bidlmensional histogram crossing two spec 
bands. An nSS4/MSS7 bidlmensional histogram is shown on 
Fig. 2.4.H in which each pair of radiance values is assigns 
to a group when a one to one relationship exists. A bispect 
representation being unable to take into account the divers 
of spectral signatures Cln four spectral bands] there are 
some pairs of radiance values which are assigned to several 
groups at the same time. Such ambiguities are pin-pointed o 
the display by a vertical bar symbol. 

Taking account of these results, it. v>ias decided to pursue ti 
subdivision of groups 1,2 and 3. 

b - Results of the subdivision phase 


The histograms of class 1, 2 and 3 were divided into 5 
equipopulation intervals. The entire process was star'ited ag. 
seperately for each group. Hereafter, we restrict our repre- 
sentation to the subdivision of group 3. 

'Correspondence analysis of the spectral signature of the pi; 
Q'F group 3 was followed by an ascending hierarchical classli 
cation whose dendogram is presented in Fig. 2.4.1. The simi- 
larities between transformed bands lead to the determinatior 
of three spectral sub-groups. 

After clustering, the individuality of each sub-group can fas 
assessed by the unimodality of its spectral histogram in iiSE 
CFlg 2.4.J). Separability between sub-groups .is however bet 
understood using the whole set of spectral bands. The averaj 
spectral signatures CFig. 2.4.K3 can be used for this purpos 
One could note that for these objets. separability is maximi 
in HSS4 and I4SS5 but very low in MSS7. This can be noticed t 
the overlapping of NSS7 histograms (Fig 2.4.J1. 

In the .display generation CFRACI) composite maps can be 
produced in which pixels assigned to each sub-group are 
represented by di-Fferent symbols. Those maps were interacti- 
vely generated on the display console of the FRACAH system. 
Fig. 2. 4. L to Fig. 2. 4. N are illustration to the subdivision 
of groups 1, 2 and 3. 

c - Interpretation phas e 

Each class now remains to be interpreted by taxonomic post- 
assistance and the objects it represents must be identified. 




Fig. 2.4.H - Joint hSS4/nsS7 histogram. Symbols are assigned to pair of radiance values 
which were associated to. specific groups in the first phase of treatment 





Fig. -2.4.1 - Ascending hierarchical tree representation 

of coded equivalent spectral bands determined 
by the subdivision of class 3 histograms 




Class 5 shows areas of clear water a indicated both by its 
spectral signature and its location. These waters 
are found to the North of the studied area : 
satellite overpass tooK place at mid-flow. The sou- 
thward tidal current through the Fosse du Fain carr 
relatively clear water across the Gois ford. Clear 
water also occurs : 

- extensively : 


. in the South part of Bourgneuf Bay where pools 
between emerged banks hold relatively calm water 
which is low in sediment content. 

. in the open sea off the Fromentlne channel where 
waters are clear at mid-flow. 

- locally : (as isolated pixels) 

. in the Trou de Sebastopol (line 94, column 19), 
a pool in a dip formed by marine erosion during 
a storn which breached the main dyke. The breach 
was repaired by rebuilding the dyke round this 
dip, outside the polder. 

. on the continent. In the pond of the new polder 
at Bouin (line 57, column 133). This pond formed 
in the lower part of this polder between the old 
seadyke which was destroyed by a storm in 1940 
and the new dyke rebuilt in 1965. 

I 

In these last two cases, the pools are so small the 
only one pixel was correctly identified as clear, 
stagnant water. Surroundings pixels were only parti 
covered in water and the classification process 
misplaced them in another class. i 

Class 4 is sediment laden sea-water 

Class 22 is associated to very turbid or shallow water. 

In this case, the refletance spectrum is influencec 
by the sea-bottom : this occurs, along Notre-Dame 
de ilonts point, in the outer tidal delta at 
Fromentine. This class also outlines the border 
between two environments. Thus, the bridge linking 
the island of Noirmoutier to the continent at 
Fromentine shows up as a row of class 22 and is 
mixed with pixels of intermediate spectral signatur 
overlaying two different environments. 

Class 21 maps the emerged tidal flats at the time of overpas 
These tidal flats line the western coast of the 
island of Noirmoutier as well as the continent. The 
also cover the top of sandbanks in the South of 
Bourgneuf Bay. Class 21 is also found in ploughed 
land or salt marsh areas physiographically similar 
to tidal flats although they fall into a different 
land use category. Finally, this class marks the 
border between classes 22 and 11 with pixels of 
intermediate spectral signature overlaying two 
different environments- 
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Class 11 shows dry area of a dominantly mineral nature : 
dry sandy beaches of Woirmoutier and the Monts 
coast ; sand or shell banks in resses parts of the 
Bourgneuf Bay coast ; sandy soils in the marsh and 
the island of Noirmoutier s sand and stone smbankmenr 
of the Bee de I’Epoids pit s and several slighty 
urban areas : Fromentine, parts of La Barre de Honrs, 
I'Epoids, la Fosse etc ... 

* I 

Glass 31 identifies vegetated areas such as the western part 
of Honts forest (maritime pines), marsh channel 
flats etc . . . 

Class 32 is similar in content and located nearby class 31 : 
Honts forest (maritime pines, holm oaks) etc ... 

Classes 12 and 13 

show cultivated fields, but the crops in class 13 
are partly bare (sandy plots on Noirmoutier Island) 
whereas the crops in class 12 have more coverage : 
on the Honts coast (near La Barre de Honts) or at 
Crosniere. 

Class B Shows grassy areas : pastures in marshes, coastal 
dunes with psammophyte cover. 

Class 33 shows marshes and damp border areas . 

Its by combining a specific ground truth experience with 
a general knowledge of spectral signatures that the operator 
can give a meaning to the classes previously defined. For this 
final step, which takes place after classification, we shall 
use the term of taxonomic post-assistance. 
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Fig. 2.4.L - Cartography of subgroups 11, 12 and 13 
obtained by subdivision of group 1 . 
Taxonomic interpretation follows : 

11 - Built up areas, beach, dry sand 

12 - Farming on calcareous, sandy soil 

(doirmoutier Island, Crosniere) 

13 - Farming along the Barre de Monts 




Fig. 2.4.M - Cartography of subgroups 21. 22 obtained by , . 

subdivision of group 2. Taxonomic interpretation is : 

21 ' Tidal flats 

22 - Very turbid water 
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Fig. 2.4.-N - Cartography of subgroups 31, 32 and 33 
obtained by subdivision of group 3. 

Taxonomic interpretation is : 

31 - Fallow marshes,- west of the Barre de Monts forest 

32 - Damp border areas 

33 - East of Barre de Monts forest 
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2. 5 - AUTOMATIC CARTOGRAPHY 


The results from all the classification techniques dsvelopped in the 
FRALIT experience can be output as colour maps an an ink plotter, by using 
the programs described .here. 

2. 5.1 - The FRACARTO program (Fig. 2. 5. A3 

This program basically maps a pre-assisted classification 
of multispectral data using spectral paramoters selected and 
input by the user.' It replaces the FRACARTE program with ' 
increased possibilities. It can also be used to map classifica- 
tions resulting from other programs such a FRACAM (c.f. 2.43. 

The main possibilities offered are ; 

“ Simultaneous processing of data from two (or more] different 
Images using all 8 Cor more] spectral bands -?or LANDSAT. 

- Rav; or smoothed data output, using a majority vote type 
decision in the matrix surrounding each pixel. ; :this process 
can be iterated, and/or used conditionally. 

- Mapping; of each pixel by a given colour and symbol or 
contouring of areas with identical taxonomic content. 

- Attainement of the best possible drafting speed according 
to the type of output, with automatic and permanent control 
of ink levels in each pen. 

The user first selects the digital image, or images, to be 
processed, and defines the scale and boundaries of his area 
of interest. 

If the proposed map exceeds the dimensions of the plotter, a 
sampling procedure is provided, enabling one in n lines (and 
one in n columns] to be mapped. 

This possibility can also be used to avoid the sixth-line 
banding often found in LANDSAT data. The map border is 
annotated with line and column tick-marks. In order to save 
time and lower costs, it is possible to leave blank an area 
of little interest within the map boundaries by specifying 
the appropriate lines and columns. 

This program is currently implemented on the CNRS IBM 370/173, 

' computer using a three-pen BENSON plotter in the offline mode. 

Pen selection is programmable and maps use up to B colours. 

For offset printing purposes, an option in the program enables 
a separation to be made of each colour overlay. These are drawn 
in black on a stable base (Mylar] to conserve maximum accuracy 
in the preparation of the offset plate. 
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Fig. 2. 5. A - rialn functions of the FRACARTO 
program,. 

Dotted paths indicate ’-possible iterations 
□f the smoothing process. 
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Whgn used for a pre-assisted classification, FRACARTO annotates 
the map with a key giving spectral characteristics of each class 
and the corresponding symbol. In every case, each BymbQl 2 is 
followed by the number of pixels, the area covered in km , and 
tv;o parallel caption lines Cwhich can be bilingual). 

Raw data mapping often shows up striping, wrong data values or 
immediate pixels with spectral signatures mixing the values of 
contrasted terrain units. The latter is a major cause of mis- 
classification, but all three effects can be avoided to some 
extent by smoothing the data before map edition. 

As the bad pixels can be differentiated from the surrounding ones,, 
smoothing requires a form of neighbourhood analysis. Optionnally, 
two types of smoothing are available : 

- maj ority-vote : here a pixel is replaced by the most 
frequent one in the surrounding neigbourhood matrix. The 
user selects the size of the neighbourhood matrix (3x3 or 
5x5) and, if required, the number of iterations. This 
techniques is fairly efficient, but eliminates much of 
the fine detail In an image (Fig. 2.5.B and 2.5.0 . 

- conditional ; this was developed to eliminate intermediate 
values. I4SS7 spectral signatures are compared for the eight 
pixels surounding the one being smoothed. Striping is 
detected by comparing the current line with the ones 
immediately before and after it. If these are identical • 
and the current pixel slightly different, its value is 
replaced by majority vote. The resulting map is more 
accurate than simple majority vote (Fig. 2.5.D) 

The results of classification smoothing can be compared to the 
corresponding unprocessed classification (Fig. 2.5.E). 

After the optimal smoothing phase, mapping Is carried out colour 
by colour, in one of three modes : 

Alphanumeric : each pixel is individually identified by an 

alphanumeric symbol. This is cos'tly in plotting 
time as each symbol requires one or more pen 
up-down calls. 

FRAMOT : a set of symbols has been selected to give varied- 

tonal values and speed of execution (Fig. 2.5.F). 
This has been found most useful both for variety 
and efficiency of use. 

Contour : identically classified pixels are oulined in groups. 

This gives the quickest results but the maps are 
only satisfactory if the classification is 
previously smoothed to eliminate isolated pixels. 
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rig. 2 . 5 .B - fcxarnplG dF majority vote smoothing with one iteration. 

nap o-F Frcmentins area at a 1/50 ODQ scale generated from 
LANDSAT2 data acquired on 29 July 1975 CE2188-101 81 1 
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Fig. 2.5.C - Example of majority vote smoothing with two iterations 

Map of Fromentir.e area at a 1/50 000 scale gensrated from 
LANDSAT2 data acquired on 29 July 1975 CE21 86-101 31 3 
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Fig. 2.5.D - Example of conditiopfil smoothing with two iterations 
map of Fromentine area at a 1/50 QOQ scale generated 
From LANDSAT2 data required on July 1975 (E218B-10181 ) 
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Fig. 2.5.E - Example of classification results before srrioothing 
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2. 5.2 - Th9 FRATHERM program 


FRATHERM is a version of FRACARTO adapted for mapping visible 
and infra red data from the NOAA VHRR. Two versions arc presently 
in use : 

The first one is a FRACARTO adapted to the output from the TSS 
(Sea Surface Temperature) programme. The two channels available 
can be combined and the graphic options available with FRACARTO 
are available. 

The second one has been perfected to make the best possible choic 
of colouror B/W scales for single-channel mapping. I/O transfers 
and computer time are reduced to a minimum [Fig. 2.5.G). 


Well designed, legible map output is the standard format for FRALIT 
results. Experience shows that other types of output (line printer 
cartography, tabular data, etc ...) are basically unattractive to potentia 
j users. On the contrary, a properly designed map. scaled for compatibility 
with standard topographic map overlays has high user value both for resear 
(geography, biology, sedimentology. oceanography) and applications (urban 
or country planning) . For these reasons a great deal of effort has been 
put into the development of sophisticated automatic napping procedures. 






ORIGINAL PAGE IS 

OF POOR QUALITY 


FRATHERn 

LIGNES 

COLOMNES 

ECHELLE: 


- UERSION 2.2 
21 A 200 n 

51 A 320 /I 

1. ’1000000 


JANUIER 1933 


rEnPEHATURE EN DEC. CELSIUS 


fig. 2.5.G - NOAA thermography generated by the fratherm program 
Satellite overpass was on OB/OQ/75 at 9H 30 MN TU. 
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2.6- GROUND TRUTH RADIOMETRIC HEASUREIIENTS 


Object reflectances measured by each of the four spectral bands of the 
fISS system build the LANDSAT spectral signature of this ODject, 

The Use of an EXOTECH radiometer (mcdel 100-A3 helps in undertaking 
ground truth operations as well as'the necessary measurement of 
atmospheric parameters. It can be used for ; 

- tne prediction of object LANDSAT signatures by "in situ’* 
measurements 

- the transformation of measured radiance by the NSS' system 
of the LANDSAT satellite into a ground albedo which is 
independent of sun angle as well as atmosphere parameters. 

2. 6.1 - The EXOTECH radiometer and its control equipment 


The EXOTECH radiometer, model 100-A, acquired by Ecole Normaie 
Superleure Ccontrat CNES 222/75) Is a multichannel ground truth 
radiometer suitable for point by point measurement. It its 
present configuration the EXOTECH radiometer is made of a 
measurement optical head made by EXOTECH and a power source 
as well as a monitoring console built by Ecole Normaie Superieure. 

a - The EXOTECH radiometer 


If is a four channel radiometer with seperate sensor optics, 
filter, detector and analog amplifier for each channel. It Is 
capable of giving four simultaneous measurements is a non- 
sequential fashion. Simultaneous measurement In the four 
channels is a compulsory feature when sensing is to be done 
from a plane or from a ship. 

The spectral characteristics of each of the four channels 
match the hSS spectral bands of the LANDSAT satellites : 

channel 1 = band NDPF4 = 0.5 to 0.6 pm, 

channel 2 = band NDPF5 = 0;B to 0.7 pm,' 

channel 3 = band NDPFB = 0.7 to 0.6 pm, 

channel 4 = band M0PF7 = 0.8 to 1.1 pm. 

The equipment is sold with 3 switchable types of sensor 
optics with the following fields of view : 


-Field of view 

Aperture shape 

Sensor optics 

1“ 

square 

lens 

15° 

circular 

diaphragm 

2ir steradlan 

hemispheric 

diffuse surface 
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The hemispheric sensor with a 2n sfceradian aperture enables 
the measurement of half-space Illumination (in watts/cm^] 
and is suitable for the measurement of down-welling solar 
illumination. The 15° circular aperture sensor enables the 
measurement of object radiance (in watts/cm^.sr) and is 
suitable for ground truth collection. The 1° square aperture 
sensor measures the radiance of a target 260 feet wl'de ’ 

(260 ft = 79.2 m3 when mounted in a plane flying at an 
altitude of 2300 maters. This , target is identical in 
geometry to a LAflDuT.T pixel. 

The output signal of each of the four channels is a' continuous 
voltage bounded by 0 and 5 volts. For each channel, four 
different gain factors can be selected. They are respectively 
of ! 1, 5, 25 and 125. 

The transformation of output volts into radiometric units 
is done by the use of a calibration factor which can take 
one of twelve values corresponding to a specific association 
between a channel and a type of sensor. Let us point out that 
the calibration o_f this type of radiometer sensing in the 
visible part of the electromagnetic spectrum is a difficult 
task. It is highly dependent on the precision of the reference 
source used during the calibration phase. 

Consequently the precision uncertainty of the radiometer is 
no better than + 5 %. The long term variations in detector 
characteristics must also be taken care of by frequent 
recalibration. 

b - Power source and monitoring console 


This console is made of a separate equipment box vsihich helds ; 
four micro-voltmeters allowing the simultaneous display of 
the four channels output, a Cadmium. Nickel accumulator power 
source and a built-in charger. A fifth commutator on the 
monitoring panel can be used in order to check on the 
accumulator charge. 

2. B.2 - Atmospheric parameters measurement 

The EXOTECH radiometer (model 100 A) is suitable for the 
measurement of atmospheric parameters such as : atmospheric 
turbidity, water concentration as well as other absorbent materials 
These measurements are -necessary in order to transform the Bateliii; 
measured radiances into ground albedo independent of solar angle 
and atmospheric parameters. This method has been introduced by 
R.H. Rogers and K., Peacock : A Machine Processing of ERTS and Ground 
Truth Data j Conference proceedings of Machine Processing of 
Remotely Sensed data [LARS] . October 1B-18th 1973 at Purdue 
University - IEEE Catalog n° 73 CHD 834-2GE, pp 4A14-4A27. 




Fig. 2. 6. A - Diffusion toward radiometer and satellite 
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If ground reflection is perfectly diffusej satellite measured radiance 
at nadir angle can be expressed as : 


L 



ir 


E T. 



( 1 ) 


where : 

L = tneasursd satellite radiance 
R_ = ground albedo 

•I 

E = irradiancB of a perfectly horizontal ground target 

= atmospheric transmittance from the target to the 
sensor per unit of air mass 

L = atmospheric path radiance 
A 

The radiance can be decomposed into two terms : 

E = E T™ cos I + E„ (2) 

□ A o D 


where. : 


E = solar irradiance outside atmosphere 
□ 

Ej^ = diffuse sky irradiance 

I = sun zenith angle 
° 1 ■ 

m = — — is the air mass along I direction 
cos I ° o 

o 

Thus, the relationship between L and is expressed by the following 
formula : 


[L 






T; (E cos I + E_3 
A 0 A 0 D 


(3) 


The different terms of relation C33 can be estimated by 'radiometric 

measurements at the time of satellite overpass. Solar angle and air mess 

m are easy to determine. Diffuse irradiance Ep from the sky is measured by 

solar occultation . As E and Ep are measured and is Known from standard 

models, T can be de' determined by equation C23. 

A 

The diffuse air radiance is estimated from the measurement of sky 
radiance Lg so that the angle between viewing direction and solar path 
Hgpj is equal to the angle between the viewing of sun and satellite Hgg. 

Diffusions towards radiometer and satellite are proportional in such a 
system geometry, the corresponding ratio being a function of the difference 
in optical path for the two directions (Fig. 2. 6. A3 : 




1 - 


T 


m 

A 



B2 


This method is only applicable for zenith angles smaller than 45°. 
When the sun' zenithal angle is larger than 45°. viewing of the sun has 
to be done for smaller than 45° and subsequently corrected. 

The main interest in such an approach Is to reduce all satellite 
measurements to physical albedo values. 5ut it is also a fact that 
in such radiometric schemes, good calibration is a critical factor. 

2. 6.3 - Ground truth measurements 


The EXbTECH radiometer Is also suitable to ground-truth measurements 
which aim at comparing satellite measured radiances and ground measured 
albedo . 


This method has been introduced by H.W. Smedes, R.L. Hulstrom and 
K.J. Ranson : The mixture problem In computer mapping of terrain : 
improved technique for establishing spectral signatures, atmospheric 
path radiance and transmittance. Houston, NASA Earth Resources Survey 
Syposium, June 1975, First comprehensive symposium on pratical applicatlor 
of earth resources survey data, volume 1-B. pp 1099-1159. 


If ground reflection is perfectly diffuse, a satellite viewing at nadir 
angle would measure tiie following radiance [equation 2} : 




and the target reflectance is expressed as : 


Kt " CL - l;) 


IT 


H S T, 


C5) 


If reflectance of a set of targets is measured at ground level or at 
a small altitude, the comparison of the values of and L shows a 
linear relationship : 


L = a X R.J. + b (B) ■ 

The intercept is equal to the atmospheric diffusion ; b = and the 
slope of this line is related to the downwelling irradiance and the 
atmospheric transmittance _ are presented on Fig. 2.6.B. 

^ ~ ir 

Consequently the spectral signature of an object in albedo related 
to the four spectral bands of the LANDSAT MSS system can be determined 
as follows : 


R. 

1 


= CL. - 

1 


Ai 


Ai 


i e i! 


where I is the set of the four MSS sensors 
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^ o-fromv H.W. SWEDES et aJ. 


Fig. 2.6.B - Relationship between measured reflectance 

at ground level and nadir radiance measured by 
satellite 
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This method has several positive aspects ; 

- the value of the measured atmospheric radiance is a 
direct contribution to the signal sensed by the 
satellite radiometer, 

~ measurement o-f atmospheric effects is done in a global 
v^^ay (downwel] ing solar irradiance and atmospheric transmiss 
effects are Integrated in the same term), 

“ no specific ground equipment set-up is necessary. 

However two type of constraints are present : 

- position of ground targets has to be Known with a 
high degree of geometric precision, 

- measurements usefulness is related to the quality of 
ground measured reflectances which implies a good 
knowledge of temporal variations of object signatures. 

2. 6.4 - Typical ground measurements 

Several "in situ” radiometric measurements were conducted in 
various- environments : 

- on the tidal flats and marine marshes of the Mont St- 
Hichel bayj 

- in harbour waters, on beaches and coastal ponds along 
the coast of Golfe du Lion, 

- in the Fromentine test area, 

- on the lawn of the Scientific campus of University of 
Villeneuve d'Ascq , near Lille, 

- on snow ccwerat the Ecole Normale Superleure in Montrouge, 

- on several samples taken from the beaches and littoral 
sand dunes of the coasts of France [measurements were 
made in natural sun light) . 

These measurements were performed in spectral bends matching 
these of the LANDSAT HSS system, using two types of radiometers 

- 820 - 10 A ERTS radiometer head manufactured fay Gairma 
Scientific, 

- EXOTECH 100 - A ERTS radiometer manufactured by EXOTECH 
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Each measurement sequence is made of 


- reflectance R. 

X 


R. = CL . . G J X 1DD / (L . 

1 Cl Cl VI 


- normalized reflectance 


"Ni ■ \ '’l 


1 « I 


with : 


G .] 
VI 


1 Cl 


= voltmeter output when measuring the target radiance 

= gain set when viewing the target 

Ly = voltneter output when measuring the downwelling 
irradiance 

Gy = gain set when viewing the sky 
i = channel index 

rieasurements are oresented on the followine tables flab. 2. 6. Cl 
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Tab. 2. 6.C - Examples of spectral signatures 

e = Laboratory sample studied in daylight 
d = depth of Secchi disk (in meters] 

•E = EXOTECH radiorneter, model 100-A 
G = Gamma Scientific radiometer, model B20-10 A 


Sample 

Haterial 

‘ Site 

MSS4 

□1 E 

beach sand 

Hont St -Michel 

31 



Bay 

1 

02 E 

beach sand 

Mont St-Michel 

52 


with shells 

Bay 

1 

03 E 

shell bar 

Mont St-Mlchel 

67 



Bay 

1 

04 E 

very dry 

Mont St-Michel 

32 



. tidal clay 

Bay 

1 

05 E 

dry tidal 

Mont St-Michel 

28 


clay 

Bay 

1 

06 E 

humid tidal 

Mont St-Michel 

24' 


clay 

Bay 

' 1 

07 E 

Salicorne 

Mont St-Michel 

.11 


- field - 

Bay 

- 1 

mam 

Obione 

Mont St-Michel 

12 

■■ 

field 

Bay- 

1 


reflectance (in, %]• . 

nssh HSS6 j nss~ 

34 42 

1.10 1.35 

55 65 

1.06 1.25 

69 64 02 

1.03 1.25 1.22 

30 37 36 

0.94 1.16 1.13 

29 34 30 

1.4 1.21 1.07 

22 . 27 

0.92 1.13 

10 53 

0.91 4.82 

10 63 

0.83 5.25 
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Tab. 2.6.C [Suite] 


Sample j 

miQH 

Site 

F 

MSS4 

Ref lectai 
MSS5 

nee tin % 
MSS6 

) 

MSS7 

09 

E 

beach sand 

Saint 

63 

65 

79 . 

61 



' with shells 

Enogat 

1__ 

1 .£3 

1.25 

1.29 

10 

eE 

i 

beach sand 

Deauduc 

42 

42 

54 

54 




point 

2 

1 

1.29 

1.23 

11 

eE 

beach sand ' 

Soulac 

64 

70 

86 

68 





1 

1.09 

1.38 

1 .38 

12 

gE- 

damped 

Soulac 

1 

38 

38 

50 

50 



beach sand 


1 

1 

1.32 

1 .32 

13 

eE 

dune sand 

Arcay spit 

56 

67 

89 

92 - 





1 

1.16 

1.53 

1.59 

14 

eE 

beach sand 

Arcay spit 

49 

56 

76 

77 





1 

1.14 

1 .55 

I 

1.57 

15 

eE 

dry ' 

Del marsh 

14 

15 

39 

50 



pasture - 


-.1 

1.07 

2.79 

3.57 

16 

eE ' 

dune sand 

Kersidan 

69 

76 

93 

93 





1 

1.1D 

1 .35 

1.35. 

17 

eE 

dune sand 

Trevignon 

49 

61 

90 

93 





1 

i:24 

1.64 

1.90 

18 

eE 

dune sand 

Rosplco 

45 

54 

82 

69 





1 

1.20 

1.B2 

1.98 

19 

gE 

dune sand 

Kerboven 

53 

61 

90 

96 





1 

1.15 

1 .70 

1 .61 
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Tab. 2. B.C (Suite) 


Sample 

Material 

Site 




m 

2D 

bE 

beach sand 

Guimorais 

55 

SB 


81 





1 

1,20 


1 .47 

21 

eE 

beach sand 

Kelcidan 

06 

BB 

81 

81 





1 

:1 

1 .23 

, 1 .23 

22 

G 

snow 

Montrouge • 

70 

78 

68 

62 





1 

1.11 

0.97 

0.89 

23 

G 

snow 

Montrouge 

58 

BS 

B4 

60 





1 

1 .14 

1 .10 

1 .03 

24 

G 

snow 

Montrouge 

60 

70 

68 






1 

1.17 

1.13 


25 

G 

thawed 

Lille 

22 

23 

68 

9 

90 



lawn 


1 

1.05 

3.09 

4.09 

26 

G 

lawn 

Lille 

9 

10 

42 

57 





1 

1.11 

4. 67 

6.33 

27 

G 

beach = 

Collioure 

14 

17 

25 

21 



pebbles 


1 

1i24 

1 .79 

1 .53 

28 

G 

beach 

Sets 

18 

15 

32 

36 



pebbles 


1 

0.83 

1.75 

1.90 

29 

G 

dry beach 

Collioure 

18 

21 

20 

20 



sand 


1 

1 .19 

1 .14 

1.13 

30 

G 

dry beach 

Port la 

22 

28 - 

30 

33 



sand 

Nouvelle 

1 

1 .30 

1 .37 

1,52 
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Sample 



32 G 


33 G 


34 G 


35 G 


36 G 


37 G 


38 G 


39 G 


40 G 


41 G 



Material 

Site 


reflectance (in % 

) 

.nss4 

MSS 5 

MSS 6 

riss7 

dry beach 

Fontignan 

15 

23 

24 

26 

sand 


1 

1 .56 

1.58 

1 .71 

dry beach 

Fontignan 

10 

13 

14 

16 

sand 


1 

1 .32 

1 .41 

1.58 

rippled dry 
beach sand 

Port la 

24 

31 

29 

36 

Nouvelle 

1 

1 ,32 

1721 

1.50 

beach sand 

Port la 

.8 

11 

11 

17 

under film 
of water 

Mouvelle 

1 

1.40 

1 .44 

2.13 

beach sand 
Under film 

Sete 

16 

24 

1.51 

26 

1.62 

27 
1 .67 

of water 


1 

beach sand 

Sets 

17 

25 

26 

31 

Under film 
of water 


1 

1 .46 

1 .61 

1 .79 

rippled damp 

Port la 

11 

12 

16 

26 

beach sand 

Nouvelle 

1 

1.05 

1.39 

2.29 

damp beach 

Port la 

16 

21 

21 

25 

sand 

Nouvelle 

1 

1 ,33 . 

1.38 

1.G3 

damp beach 

Sete 

25 

34 

37 

38 

sand 


1 

1 .32 

1 ,46 

1 .50 

damp beach 

Frontignan 

13 

25 

25 

30 

sand 


1 

1 .33 

1 .33 

1 ,61 

shell sand 

Sete 

27 

37 

42 

50 



1 

1.39 

1.57 

1.87 

shells and 

Sete 

36 

45 

47 

8B 

pebbles 


1 

1 .26 

1.32' 

1 .86 













































































7D 


Tab. 2. S.C [Suits) 


Sample 

Hatsrial 

Site 


mam 



43 


shells 

Sete 

24 

39 

45 

52 





1 

■ 1.64 

1.89 

2. 18 

44 

G 

. 

shells 

Frontlgnan 

14 

20 

24 

27 





1 

1.40 

1 .67 

1.89 

45 

G 


Sete 

10 

14 

17 

21 



sand 


1 

1 .34 

1 .67 

1 .99 

46 


muddy 

Sete 

12 

18 

19 

21 



sand 


1 

1 .43 

1 .57 

1.73 

47 

G 

muddy 

Sete 

13 

18 

19 

24 



sand 


1 

1 .39 

1.48 

1.63 

48 

G 

muddy 

Sets- 

15 

21 

21 

24 



sand 


1 

1 .44 

1 .45 

1 .59 

49 

G 

muddy sand 

Sets 

18 

23 

26 

31 



with film 
of water 


1 

1 .27 

1 .47 

1.72 

50 ' 

G 

red soil 

Frontignan 

10 

30 

43 

47 





1 

2.94 

4.18 

4.60 

51 

G 

muddy sand 

Sete 

7 

10 

14 

10 



with film 
of water 


1 

1 .33 

1.94 

1 .45 

52 

G 

muddy sand 

Sete 

7 

10 

14 

10 



. with film 
of water 


1 

1 .25 

1.86 

1.35 

53 

G 

muddy sand 

Sete 

8 


15 

12 



with film 
of vjater 


1 

mm 

2.01 

1 .54 
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Material 


□bione -field Sete 


□blone field Sete 


salicorne 


field 



reflectance (in 
MSS1 f MSS5 i MSS6 { MSS7 


12 41 59 


1.50 5.31 7.53 


I2 42 BD 

1".72 5.84 8,38 


5 12 2D 

1.1S 2.86 4.93 


d = 3,0 m 


Port V/endres 1.86 0.81 0.41 0.27 

1 0.44 0.22 0.15 


d = 3.0 m 


d = 3.5 m 


ocean 
d = 3.5 m 


ocean 
d = 4.0 m 


d = 4.5 m 


Port Vendres 


Port Vendres 


Port Vendres 


Port Vendres 



0.45 0.30 

0.24 0.18 


0.30 0.21 

0.18 0.12 


Port Vendres 2.15 1.13 -0.45 

1 0.53 0.21 




0.3B 0.21 

0.19 0.11 


0.44 0.22 

0.23 0.11 


d = 4.5 m 


Port Vendres 2.13 1.25 0.75 0.85 

1 0.59 0.35 0.40 


ocean 
d = 3 .5 in 


d = 4.5 m 


Port Vendres 


Port Vendres 



0.28 0.10 
0.15 0.10 
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2. 7 - DETERniNATIDIN OF TIDAL CURRENTS AND TIDAL HEIGHTS 


Since a large part of the works of the FRALIT team are related to the 
tidal coast of France, it was necessary to compote the tidal heights 
as well as the speed and heading of the tidal currents at satellite 
overpass. Those determinations were drawn from data collected by the 
Service Hydrographique and Oceanographique de la Narine and sometimes fi 
other data sources CE.D.F., direct observations, 

Figures 2.7.1 through 2.7.13 provide this tyne cf information for 
15 images of the French Atlantic littoral which were computer processed 
processed by the FRALIT team ; those images are ordered by geographical 
area and for each geographical area by chronological order. 
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3 - SIGNIFICANT RESULTS AND ACCOnPLISHHENTS 


3. 1 - AN AUTOMATIC CLASSIFICATION OF LANDSAT-2 DIACHRONIC DATA ON THE 
FORESHORE OF THE ISLAND OF JERSEY 


3. 1.1 - Data description 


Landsat-2 data from June Sth 1975 and July 30th 1975 were used 
to drawi up a diachronic [multitemporall map of the island of • 
Jersey. These two dates were chosen for the differences in tide 
levels between them. In the first case, satellite overpass, took 
place 20 minutes after low water at St-Helier, and tvjenty 
minutes after high water in the second case. 

3. 1.2 - Data processing 

Diachronic processing requires the overlayging of the two 
digital images, using significant landmarks identified on 
single channel imagery as tiepoints. At the end of this step, 
a single date file is available, , where each pixel is characterlar 
zed by eight different spectral bands [4 hands for each date] . 

Preliminary to dat analysis, radiance histograms in each of 
the B spectral bands are divided into radiance intervals. This 
subdivision Is used in order to generate a binary coding of the 
original data. When the radiance histogram of a spectral band is 
unimodal it is divided in isopopulatien intervals. However if 
, this distribution in bi or multimodal, a rapid interpretation 
of the main spectral contrasts allows to interactively adapt the 
radiance intervals to this situation. This second situation V'fes 
the case for the NSS7 bands at the two differents dates analysed 
here. The emphasis was put on the contrast existing between sea 
and land signatures. 

After correspondence analysis of coded data using the FRACAM 
system, a hierarchical classification reveals the similarities 
existing between spectral bands. The groups of spectral bands 
thus determined are displayed as nodes of a hierarchical tree 
structure (dendogram) in Fig. 3.1. A’. 

The lenght of each branch of this tree is proportional to the 
distance between the groups it links. Operator interaction is 
required for the determination of the maximum group diameter 
as well as the number of desired groups. This step is followed 
by an adaptative clustering of radiance spectra leading tc a 
spectral classification. 

It Is possible to control the results of this classification fay 
displaying the class spectral histograms and spectral statistics. 
The analysis can then be refined in a layered fashion by reprocss 
sing each class for further subdivision. The present study being 
mainly related to the tidal flats dynamics, classes 3 and 7 were 
subdivided in 3 subgroups each. 
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3. 1.3 - Interprstation of classification resul-ts 


The map legend CFig. 3.1.B3 shows each taxon In its relative 
■-■position within the hierarchy defined by the classification. 

, I ^ 

For interpretation purposes, two types of taxons have been defined 
according to their variability in time. On the whole, taxons with 
a similar spectral signature at both dates can be opposed to those 
with strongly varying spectral signature. 

a - Categories with low diachronic variation 


. Class 1 shows the upper part of beaches [South and south-east 
coast) v/hich was not covered by the neap tide of July 30 th. • 
It also borders class 2 and occurs inland, generally, vjhere 
it shows low-density built-up areas. 

, Class 2 covers the town of St-Helier and the discontinuous 
urban areas along the South and East coasts. 

. Class 31 occurs along the North coast in moorland. This subset 
of the vegetation class was individualized during the second 
stage of classification. 

. Class 33 shows low density built-up areas. 

, Class 4, S and 6 are different types of vegetation. The same 
symbol has been used to map them because of the lack of 
approriate ground-truth data for accurate interpretation. 

. Class 71 represents bare rock, in the cliffs to the North East 
of tlie island, and in the new harbour area under construction 
southward of St-Helier. 

. Finally class 8 is an example of "false stability" in time 
since it shows areas which were covered by the sea at both 
dates . 

b “ Categories with strong' diachronic variation 

These are mostly foreshore areas, with a land-type spectral 
signature at low tide and a sea-type signature at high tide. 

At the first level of classification, the foreshore was 
defined by class 3 and 7, but some of their subclasses are 
in fact land area as mentioned above. 

Class 31 and 72 show different moisture levels ideally 
illustrated in St-Aubin's bay where class 1 , 31 and 72 occur 
in sequence between the high and low tide marks. 

Class 71 covers the emerged parts of the rocky shoals south- 
east of the island. 
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Fig. 3.1.B - Automatic diachronic classification of Jersey 

island using LANDSAT-2 




3. 1.4 - Critics 


The diachronic classification used for this study has both the 
simplicity and rapidity of implementation of post assisted 
automatic classification, and the good identification obtained 
by diachronic techniques. For observation of rapidly changing 
terrain types such as the foreshore, diachronic techniques make 
the- best possible use of the respetirivity and fairly good 
geometric repeatibility of satellite remotely sensed data. 
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3. 2 - INTERPRETATION OF LANDSAT DATA OM THE LANDWARD PART OF THE 
FROMENTINE TEST AREA 


3. 2.1 - Data description 

Data analysed on the Fromentine test area- are related to the 
following dates of overpass : 

27 September 1972 - LANDSAT 1 - 

11 July 1973 - LANDSAT 1 - 

29 July 1975 - LANDSAT 2 - 

3. 2.2 - Data processing 


The Fromentine test area has been analysed using various 
analysis techniques described in part 2 (Cf 2-2 to 2-5]. 

The results presented were respectively related to the 
follovjing dates : 27 September 1972 (cf. 2-3], 11'July 1973 
(cf 2-4] and 29 July 1975 (cf 2-2 and 2-5]. 


3. 2.3 - Interpretation of results before classification smoothing 


The interpretation is restricted to the landward part of the 
Fromentine area including the sand beaches which are often 
difficult to spectrally separate from the barren coastal dunes 
in the southern part of Noirmoutier island as well as along 
the Breton marsh. 

Comparisons are made between results more than between process! 
methods because they usually depend upon the quality of operate 
assistance. This parameter has obviously improved in time with 
the Knowledge acquired by the FRALIT team on the test area. 

A certain number of conclusions can be drawn : 



a - All methods of classification easily separate water from 
land categories 

b - All methods of classification map relatively abondant sand 
formations which present high radiances for the whole set 
of NSS spectral bands. Their spectral separation is usually 
the best (cf. 2-3, AGREG] . This category is associated to 
barren sands, sandy bars and various open spaces with varyl 
proportion of buildings and sands. 

c The various classification techniques are usually successfu 
in the separation of the maritime pine forest (Monts forest 

. The FRACARTE classification (Fig. 2.2.C, 29 July 1975] maps, 
in category 8 the well outlined forest area as well as 
polders which were taken over by the sea. 

. The AGREG classification (Fig. 2.3.D, 27 September 1372] 
displays a restricted and homogeneous forested area (categor 
Certain parts of the forest under peculiar stress condition 
have been omitted but on the average the recognition of the: 
Monts forest by this classification is excellent. j 
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. ThG FRACAM classification C'Fig. 2.4.Nj 11 July 19753 splits 
the forest into two parts j 

31 to which are associated the low parts of the old 

■Grande Prise poldsr and a part of the flonts forest > 

and 33 where an other part of the Monts forest is displayed 

A FRACAM classification of another date (29 July 19753 is 
presented' in Fig, 3, 2, A. It also points out a division of the 
Monts forest into two parts-. The first symbol is related to 
a dense stand of maritime pines which comprises numerous young 
trees along the Monts coastline or older patches of wealthy trees 
It is this denser part of the forest which is associated to 
the old Grande Prise polder taken over by the sea. This confu- 
sion is due to the fact that the spectral signatures of those 
two types of objects are quasi Identical [cf. 2.2.B). 

Could i-t be that the use of a multitemcoral approach v/ould 
eliminate this ambiguity which is presenr for the LANDSAT 
scenes of 11 July 1973 and 29 July 1975 ? This similarity 
between spectral signatures of otherwise distinct types of 
objects makes their possible separation a difficult task. In 
particular when this confusion 'Can be repeated in the data of 
distinct dates (11 July 1973 and 25 July 19753. The use of 
diachronic data could save this problem if for instance one 
of the survey dates would' correspond to the flood of a very 
high tide (which is impossible due to the time periodicity 
of LANDSAT coverageO or the period following the uncovering 
□f old polder by the sea. In this second type of situation the 
presence of water in one of the environmeiits Cold polder) 
associated to its absence in the other on^ (forest3 would 
help in order to perform a successful separation. 

d. - The other types of land use are mapped with a varying amount 
of success by the various classifications. We will only compare 
the results obtained for the 29 July 1975. Jby FRACARTE (cf 2.2. C3 
and FRACAM (Fig 3. 2. AO. '-o. 

’ • , * >r 

□neremark has to be done before any interpret4'tipiT..-,The mapping 
□f agriculture areas was not undertaken in full 'oetaill , during 
preassistance of FRACARTE. Motivations of the operator were 
biased by a coastal mapping approach so that more details were 
sought in tidal flats or ocean waters than in the vegetated 
marshes. On the contrary the operator while performing the 
FRACAM coassistance phase is taking .into 'account the spectral 
■ statistics of the data more than, th'e ”a priori" nature of the 
environments to be m^ped.- \ \ 

For the FRACARTE' classification two symbols'^nly'’are associated 
to the marsh land use type. ' ' 

- symbol 9 is related to psammophyte dune cover, most 

of the marsh pastures as well as some of the cultivated 
fields of Noirmoutier island and the Breton marsh j 

- symbol 10 represents essentially cultivated areas as well 
as some grass pastures. 
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On the contrary, the FRACAM classification results propose 
six different categories on the marsh : 

- symbol 3 : cultivated areas (mostly on the rim of the 

Monts forest) ; 

- symbol 4 : partially built open spaces ( Fromentlne, th: 

Sarre de Monts pier. La Fosse) j 

- symbol 6 ; psammophyte covered dunes, marsh pastures ; 

- symbol 7 : cultures with small plant canopy cover most 

on sandy soils ; 

- symbol 8 : vegetation along drainage ditches and canal 

- symbol 9 : lower parts of polders taken over by the se 

3. 2.4 - Interpretation of results after classification smoothing 

We will compare the mapping results of the raw and smoothed 
FRACAM classification of the 29 July 1975 data set. Cartograph 
improvements due to classification smoothing will be assessed 
(Table - 3.2.B) . 

a - Common properties to the two types of smoothing 

It is noticeable that conditional and majority vote smooth 
algorithms have the common property to suppress artificial 
categories generated by the classification processes at th 
boundary between two types of well contrasted objects. 
Symbol 5, for example, is mostly representative of muddy 
tidal flats but it has been introduced as an artefact alon 
the seaward part of the Monts forest sandy coast (Fig 3.2. 
This mlsclassification can be eliminated by each of the 
smoothing algorithms working on "clean" categories such as 
the sandy beach (12) and the turbid part of the sea water 

In the Nolrmoutier island, the Trou de Sebastopol (lines 
24 to 25, columns 5 to 6 in Fig 3. 2. A) is partially stable 
under a two iterations conditional smoothing. The pixels o 
this water pool take very low values in M337 which are ver 
different from those taken by the surrounding pixels. This 
spectral discontinuity insures the stability of the Trou d 
Sebastopol pixels through conditional smoothing because th 
process is more influenced by spectral contrasts than by t 
relative geographic position of the pixels (Fig 2.2.E). Th 
pool can still be outlined after a simple majority vote 
smoothing (Fig 2.2,0 but if this type of algorithm is 
iterated, the Trou de Sebastopol can no longer be identifi 
because it is merged with the nearby ocean surface (Fig 2. 
This confusion is due to the suppression of the polder dit 
which separates the pool from the Bourgneuf Bay . Majority v 
rule does in fact drastically alter linear shaped objects 
like ditches, canals, piers, bridges ... 
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Cultivated plots reprsssnted by two symbol 3 (ligne 99 and 
columns 111 to 112, for example] are surrounded by symbol 4 
pixels in Fig 3. 2. A Cunsmoothed classification]. 

Those plots are left untouched by a two iterations conditional 
smoothing due to the feet that their spectral individuality 
is good and their shape is not too similar to horizontal 
striping artefacts (Fig 2.2,E], On the contrary, the presence 
of category B on the lines 98 and 10G Pauses the removal of. 
those isolated plots at the first iteration of the majority 
vote smoothing process CFig 2.2,C], 

b “ Differences between properties of the two smoothing algorithms 


Conditional smoothing results are usually different from those 
of the majority-vote algorithm, even If both of them succeed 
in suppressing most of the ’‘artificial” category assignments. 

The numerous isolated singularities which are present in our 
landscapes, often influenced by human activities, are usually 
unaltered by conditional smoothing. This is the case of small 
entities such as : the Trou de Sebastopol pool, the small 
clear cuts cf the Ments forest, the cultivated plots scattered 
in the pastured marsh and associated to farm buildings, 
pastures intermixed with the fields along the' rim of the 
Monts forest, small pools of the active or abandoned salt 
mills. 

The majority vote smoothing, on the contrary, has a trend 
tovjarda strong generalizations. It is more suited to the 
enhancement of global contrast such as the one which separate 
the forest into two parts. The specific land use pattern along 
the forest rim, mixing cultivated field ans houses, is also 
well outlined by majority vote (Fig 2.2.D]. 

Each of the two smoothing algorithms has particular advantages 
and drawbacks, they should be selected in relation to the 
type of geographical theme to be studied. 
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Fig. 3. 2. A 


FRACAn Classification display of the Fromentino 
area for the 29 July 1975 
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FRACAN classification (29 July 19753 



ho Smoothing 

With FRACARTO smoothing 

symbol 

n° 

FRACARTO 

mapping 

program 

ilajority vote j 

Conditional 

simple 

iterated 

Fig. 3. 2. A 

Fig. 2.2.C 

Fig. 2.2.D 

Fig. 2.2.E 

1 

1 095 

1 236 

1 303 

1 124 

2 

433 

298 

256 

390 

3 

614 

472 

44B 

t'- 

541 7 

4 

472 

205 . 

146 

330 

5 

1 251 

1 260 

1 264 

1 247 

6 

3 326 

4 B12 

4 864 

4 068 

7 

967 

733 

665 

B79 

8 

1 036 

651 

501 

875 

9 

236 

163 

145 

165 

10 

3 330 

3 180 

3 200 

3 170 

11 

2 939 

2 920 

2 323 

2 341 

12 

497 

554 

569 

554 

Total 

16 800 

16 284 

16 284 

IB 284 


Table 3.2.B - Number of pixels pertaining to each category 

When smoothing data, a rim of pixels is lost 
along each sides of the processed area 
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3, -3 - INTERPRETATION OF LANDSAT DATA COVERING THE LOIRE ESTUARY 


An accurate Knowledge of the position of sandbanks in the Loire estuary 
CFig. 3.3.A1 Is essential for navigation up to the major harbour at 
Nantes - St Nazaire and the new docks at Donges f10 Km upriver] . Navigablt 
channels are frequently dredged to maintain a constand depth. Outside the- 
channels the shallows silt up 'continuously . This is most apparent off 
Paimboeuf where low tides uncover numerous sandbanks, visible on charts 
and survey maps. Extension and location of these banks varies considerably 
according to map type and date. There are three main causes for these 
variations : 



- The banks have very weak slopes and a 10 cm error in soundings or 
levelling can lead to tens of meters in planimetric error. 

“ The nature of these banks [mud or sand] makes accurate bathymetrit 
or topographic measurement practically impossible. 

- The combined action of currents, waves and estuaire sediment 
transport is constantly modifying the size and shape of each bank 
These actions can be very rapid, especially during a storm. 

The long and costly depth-sounding missions are not very frequent. Survey: 
were carried out by the Port of Nantes - St Nazaire harbour authority in 
1969 [resulting in a 1/25 000 scale map] and in 1976-77. 


Aerial surveys are also fairly infrequent, and are merely used for 
updating terrestrial topographical maps. 


knowledge of sandbank movements is thus -very inaccurate and is mostly 
based on word-of-mouth reports from local fishermen who describe the 
difficulties encountered in their navigation and the new channels they 
OQcasionnaly find. 


3. 3.1 - Landsat imagery used in this study 

Technical conditions seemed, to be ideally suited for an 
application of satellite imagery to see if this was adaptated to 
sandbank monitoring. 


The following table shows the images used and the corresponding 
tidal height. 


Figure number 

Data 

Source 

Tide height 
[meters] 

Cycle 

Fig. 3.3.D 

Sept. 27 1972 

NASA 

2.37 . 

EBB 

E 

Narch 8 1973 

It 

• 1 .44 

EBB 

F 

July 29 1975 

II 

3.18 

EBB 

G 

March 28 1976 

relespa- 

zio 

1.46 

FLOW 

H 

June 25 1976 

If 

1 .55 

FLOW 


Two other images were available covering this area but were not 
used because of their high tide level (over 4 meters] . 
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3. 3.2 “ Data processing in relation to local conditions of Interpretatlc 

Interpretation of the maps shown here is not straightforward 

a - Three of the images show an ebb tide . This situation is 
characterized by powerful seaward currents with a heavy , 
sediment load which is largely carried out to sea. These 
currents sweep the shallows and raise most of their fine 
- sediment cover. It is therefore fairly difficult to distingu: 
the actual sandbanks from tui'bid water. 

The flow tide which predominates on the remaining two images 
brings clear oceanic water into the estuary and the banks ai 
then clearly distinguishable. 

Four transects (Fig. 3.3.B and 3.3.C) have been drawn up 
near the Banc de Bilho : two show ebb tides with Nasa-type 
digital values (5 to 8 in MSS7 for emerged banks) and two 
flow tides with Telespazio digital values (20 to 40 MSS7 foi 
emerged banks). In both situationsj the distinction between 
banks and water Is easiest when the river outflow is lowest. 
During the low discharge periods, the estuarine turbidity 
area does not extend this far down channel whereas in flood 
periods this turbid area extends further seaward. 

b - The part of the river between the port of St Nazaire and thi 
new docks at Donges is heavily industrialised. Aerial pollu; 
is occasionally blovjn over the study area (cf Fig. 3.3.F) - 
has a similar spectral signature to sandbanks. This can leac 
to certain misclass'ificaticn. Urban areas are also misclass: 
fled but their known position does not hinder interpretation 

c - Finally, the most difficult element of misclassification is 
the quasi-identical signature of turbid water and sandbanks 
Fig. 3.3.E shows a long turbidity plume near Mindin Point, 
developped in the same direction as the current. This is 
probably connected with the drilling operations in progress 
then for thS' construction of a bridge across the estuary at 
St Nazaire. There is a second, lighter, plume to the north, 
where drilling was less advanced. The muddy waters move 
alternately up and downstream with the tides. 

3. 3.3 - Significant results on ‘the evolution of the sandbanks 

Observed variations cannot be accurately measured due to the . 
differences in sea level between satellite overpass (1.74 to 0.! 
meters) ! pixel size Is not compatible with this type of scene | 
variation. 

a - Overall evolution from 1972 to 1376 : Throughout this perio 
at similar tide levels, the banks reduce in area. The two 
high river discharge Images (3.3.E and 3.3.G) show over a 
two year period an identical outline for the Bilho bank to 
seaward whereas upstream the bank, has receeded in same tim 
to a line joining Paimboeuf to Montoir. The Brillantes bank 
has receeded at both ends, partly because of dredging 
operations in the access channel to Donges harbour. 
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The Nindin bank, clearly visible in 1S73, has almost 
disppeared in 1975 where it is obscured by smoke from the 
St Nazairs area. The St Brevin bank, extensive and separate 
from the coast in 1973, is not apparent in -1976. Finally the 
Pointeau bank with its typical y-shape is considerably reducac 
in 1375. Thus the overall picture is an erosion of all the 
banks between 1973 and 1976. 



The three low river discharge images confirm this interpreta- 
tion especially- since the water level is lowest at the time 
of the most recent image ; Bilho and Briliantes banks are 
ti’uncated to the north, Pointeau bank is thinned out and 
St Brevin hank has disappeared. 

b - Seasonal variations 

A seasonal cycle can be distinguished between the March 
images on one hand and the June and July ones on the other. 

. The Bilho bank is seasonally stable, with little or no 
variation in its outline but the highest , points of this 
bank [darkest tore on the maps) move fairly regularly upstrea 
from March to September. Remarkably, one of the culminating 
points in the March imagery [Fig 3.3.E and 3.3.G), located 
farthest downstream, disappears completely in the other Imagj 

The St Brevin bank Is well separated from the coast in 1S72 
[Fig 3.3.D), and only linked to the coast at Mindin Point. 
During the winter, the whole bank moves nearer the shore 
and increases in height at both ends [Fig. 3.3.E). In March 
1976 [Fig. 3.3.G) the remaining oart of the bank is up 
against the coastline for most of its length. In June 1976 
[Fig. 3.3.H) the bank is practically indistinguishable .from 
the coastline. It appears, thus, that this bank is generally 
further offshore in the springtime than in the summer. 

. The Pointeau bank forks into branches in the spring (3 in 
1973 ! Fig 3.3.E, 2 in 1976 : Fig 3.3.G} beginning at the 
Pointeau point and apparently modelled by currents. In June 
these branches are assembled into a triangular shape and 
in September a single curved ridge remains, probably reworks: 
by the swells which strike it head-on. 


c - Possible causes for these variations . For the three banks 
which show a seasonal evolution [Bilho, St Brevin, Pointeau) 
the spring state is conditionnad by the geomorphological 
action of ebb-tide currents. The banks are stretched downstreer 
into divergent branches in the outer estuary. In the inner 
estuary, the banks move downstream where they are built up to 
a peak. Conversely in the summertime, when river discharge is 
reduced, wave influence is predominant and the banks are movad 
inshore and their outline is simplified, occasionally concave 
seav;ard. The banks of the inner estuary move back upstream 
and their situation is stabilised until the next spring floods 

3. 3.4 - Conclusion 

-Long term evolution of the Loire river sandbanks has been evidence: 
despite the scarcity of available imagery. Possible explanations 
have been proposed for the evolution but they can only be confirms: 
by use of more frequent data covering other periods than spring cr 
summertime. 
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Fig. 3.3.G - Cartography of sandbanks from LANDSAT 

[March 18 19761 
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Fig. 3.3. H. - Cartography of ■sandbanRFTrQnfT-ANDSAT 

CJune 26 19763 
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3.4- INTERPRETATION OF GDTA SURVERYING OF MONT SAINT-NICHEL BAY 


3. 4.1 - DATA SET SPECIFICATIONS 


The aerial survey was carried out tay the GDTA CGroupement pour le 
Develeppement de la Teledetection Aeorspatiale) on 26 June 1975 using 
a plane of Intitut Geograohique National . 

a - Description of radiometers 


Two DAEDALUS multispectral scanners were used in the visible and the 
'thermal infra-red part of the electromagnetic spectrum. 

These electro optical sensors are equiped with a mirror rotating at 
80 r;-q.s. and the entire optical system has an instantaneous field 
of view of 2.5 milliradians . 

In the visible spectrum the DABDALUS radiometer is equipped with a 
set of ten detectors whose waveleng-bh characteristics are listed in. 
Tab. 3. 4. A. In the thermal infra-red a second DAEDALUS radiometer was 
used with two spectral bands whose sensing windows are : 

- 3 to 5 micrometers 

- 6 to 14 micrometers 

The analog output signals were recorded on a SANGAMO SABRE III 
14-channel recorder along with a digital time bass and continuous 
- flight parameters such as : heading, roll and pitch. 

b - Calibration schemes 

The DAEDALUS radiometers are equipped with two internal calibration 
sources which are viewed at each mirror rotation before and after earth 
scan, they are called positive and negative reference values, respecti- 
vely R and R~. 


Transformation from output voltages to calibrated values is done on 
the basis of a running average of the R and R~ time series with a time 
lag of 9 scans. Finally each of the 512 sample points of a scan line 
are transformed into calibrated values, using the following formula ; 


where : 


V . = (C. X - R J / (R . - R .] 
Cl 1 bi ai ai ai 


= calibrated radiometric values of spectral band i 

“ output voltage of spectral band i before transformatio 

C^ = amplitude coefficient of spectral band i used in 
order to generate an 6 bits C256 level) coding of 

R .= running average of negative reference values for 
spectral band ± 


1 + 

ai= running average of positive reference values for 
spectral band i 
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Channel 

Spectral band (in pm) 

Me De .(I) 

1 

0,38 - 0.42 

3,0 % 

2 

0,42 - 0,45 

1,3 % 

.3 

0,45 - 0,50 

0,4 % 

4 

• 0,50 - 0,65 

0,3 % 

■ 5 

0,55 - 0,60 

0,2 % 

6 

0,B0 - 0,65 

0,2 % 

7 

0,65 - 0,70 

0,2 % 

8 

0,70 - 0,80 

- 0,2 % 

9 

0,80 - 0,89 

0,2 % 

10 

0,92 - 1,10 

0,3 % 


Table -3.4. A “ Spectral bands of the visible 
DAEDALUS radiometer 


[I] : Ng Dg is the minimum detectable 
reflectance increment 
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Ths amplitude coefficient C. is determined for each visible spectral 
band i as ; 

C. = L. X K. (visible and near infra-red) 

1 11 

with 

L. = calibration coefficient efiual to the equivalent 
radiance of the built-in reference source 

= conversion coefficient ensuring an 8 bit outpuc cods 

1 

It is worth noting that K.. is different from one spectral band to the 
next but is also dependent on the surveyed scene. 

In the thermal infra-red, determination of is slightly different : 

C. = (T.„ - T.,) xK. [thermal infra-red) 

1 i2 i1 1 

with ; 

= equivalent black, body temperature of the negative 
reference source 

T ^2 ” equivalent black body temperature of the positive 
reference source 

c - Survey.' planning 


The surveyed area (Fig. 3. 4. B) is organized as a meridian transect going 
from the Dol marsh in the south to the Chausey islands in the north. 

It displays transitions between various environments such' as : a precem- 
brian formation, a cultivated marsh, some semidentary and rocky tidal 
flats and various gradations of water turbidity. 

Data gathered on 26 June 1975 are organized in five different overpasses 
labeled a to e(Tab.3.4. C) . Results presented here are related to 
overpass a because it was the only one for which photographic documents, 
DAEDALUS analog displays and corresponding digital data were available. 
Furthermore at time of survey (15h 20 mn) hydraulic conditions were good 
because tidal flats were uncovered by a low tide of coefficient D.79. 

During the days preceeding the survey no rain precipitation was noticed 
and winds were measured at the Pointe du Grouin semaphore as N.-W with 
a 3 to 4 m/s speed. 

3 .4.2 - Data processing 

" Computer compatible tapes were generated by GDTA for overpass 1 and 3 which 
are labeled a and e in our nomenclature. Data analysis was limited to a 
subset of overpass a which regrouped a representative sample of objects 
to be mapped. 

Different data processing methods have been used regarding the thermal infra 
red temperature imagery and the measured radiances in the visible and near 
■ infra-red. 
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overpass 

■ 

Hgading 

Altitude 

Photography 

Cl] 

DEAOALUS 

analog 

imagery 

DAEDALUS 

digital 

data 

a 

Cancale 

10° 

2510 m 

Yes 

Yes 

Yes 

B [4) 

Cancale 

210° 

'4875 m 




r C5) 

Mont 

St -Michel 

35° 

5060 m 




6 

Cancale 

10° 

2510 m 

Yes 

Yes 

No 

e 

CsncalB 

CD 

O 

o 

2510 rri 

No 

Yes 

Yes 


Table 3.4.C - Hissions characteristics of 

Mont Saint-Michel GDTA survey 


(1 ] t Infra-red photography taken with a WILD RC8 camera 
of 152 mm focal length and 74° aperture 

(2!) : channels 3, 4, 5, 10 and thermal T.R. 

(3] : channels 1 to 5,7, 9, 10 and thermal I.R. 

(4] : documents not available 

(5] : -documents available at l.G.’N. 
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a - Thermal infra-red processing 


Data recorded for the 8-14 ym band are transformed into surface temperc 
tures using calibration parameters stored in the first records of the 
magnetic tape. Output of this preprocessing phase in a file of pixel 
values in, Celsius degrees uncorrected from atmospheric absorption. 
Format of this file has been made compatible with the input format 
of the NCAA VHRR thermography processing program CTSS system] . 

Following treatment done by TSS is in two phases ; 

- removing of striping effect noise by resetting the sequence of 
line histograms 

“ smoothing using a 7 x 7 boxcar filter conditioned to contrast 
'enhance significant boundaries. 

Results can be displayed on various peripheral equipments : line prlnte 
electrostatic printer-plotter (VERSATEC) for intermediate maps and 
incremental color pen plotter CBENSON) for final edition. 

b - Visible and near infra-red processing 

Measurements recorded in the eight channels (Cl to C5, C7, C9 and C10] 
are preprocessed in order to estimate and subsequently correct for 
along seanlines systematic bias trend. This trend is due to the varlati 
along the scan of the relative geometry of viewing angle versus solar 
angle. 

After bias correction, the different typesof ground objects are cate- 
gorized using the FRACAM interactive classification system (cf 2.4.) 
Some groups can be subdivided after the initial classification phase 
'using the adaptative descending classification approach. 


Final results are displayed as color maps using the FRACARTO program 
;fcf 2.5). 

• 3 . 4.3 - Significant results 

a -. Test area presentation 


The. test area presented here includes most of the infra-tidal and 
continental categories which can be characterized-isln the west part of 
Mont Saint-Michel bay. 

Category and terrain description are drawn from observations which 
were made at airplane overpass time or a few hours after it. 

a - 1 - Infra-tidal zone 

The infra-tida.l zone is characterized by flat expanses always 
carbonate rich with a N.W oriented slope ranging from 2.5 0/00 
to 3 0/00 along tidal channels. 

a - 2 - Infra-tidal zone 

In front of Vilde - la - Marine village where ground- truth was 
collected, the tidal flats are wide (of the order of 3 km in 
width) and dip to the sea v^lith a gentle slope (of the order cf 
3 0 / 00 ) . Sediments are made of "tangue” which is granulometry 
and minerally defined as follows : 
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Fig 3.4.D - Localisation of categories (taxons] 

present on the infra and inter-tidal zones 
A to F : diatomee samples j 1 to 6 : sediment 
samples ; 

I : turbid to non turbid water ; 

II ; flooded mussel nurseries 

III : dry mussel nurseries ; 

IV : mussel nurseries on diatomee covered "tangue" j 

V : diatomee covered "tangue” ; 

VI : bare "tangue" ; 

VII : fish traps ; 

VIII : channel ; 

IX : calcareous sand bar ; 

X : tidal flat Cschorre] 


- particle size median between B and 103U m ; 


- percentage of the 2 to 200 Urn fraction ranging between 
51 % and 84 % 


- percentage in carbonates ranging from 26 % to 65 ^ 

Thee is a trend to larger particle sizes and coarser sediments 
near the coast line CFig.3.4.D and Tab3.4.E. samples 1 and 3). 

During the survey, concentration in organic matter was found to 
constant over the whole tidal flat surface with values ranging 
from 5 % to 7 %. On the contrary sediment water content was 
decreasing landward from the low tide line [36.5 %] to the salt 
marshes (18.5 %) . 

Finally, large expanses of the low flats were covered with a fir 
layer of dlatomees [Tab 3.4.F). 

In the upper part of the inter-tidal zone the monotony of "tangc 
flats is disrupted by sandy levees. The shell-rich mate-rial of 
these levees is quite different from "tangue" granulometry as 
well as their concentrations in water, carbonates and organic 
matter [Tab3.4.E, samples 4 and 5). 

The only main tidal channel is the Canal des Allemands. It is 
deeply cut into the upper part of the flat resolving to a modest 
opening near the sea line. 


It is also noticeable that the "tangue” flat bears a large numbe 
of human artefacts ; fish traps which are one to two meter high 
wooden constructions, mussel shells poles scattered seaward of 
the traps. 

Near the upper part of the flat, the salt marsh [schorre) is 
gradually taking over the "tangue". The salt marsh is limited ir 
the test area to a narrow triangular surface on the east bank ol 
the Canal des Allemands. Its gramelometry and mineralogy are sin 
to those of the "tangue” [Tab3.4.E, sample 6). It is, however, 
different from the barren flat in its morphology and bears salt 
tolerant plant grouth mostly composed of Puccinella maritima . 
Aster trlpolium and a thin strip of Obione po'rtulacoldes along ^ 
channel banks. ' 

a - 3 Continental zone 

Behind the Duchesse Anne Lyke, the test area extends over the i 
Dol marsh. This cultivated area has soil properties which are | 
very close from "tangue". 

Results Interpretation 

b - 1 Thermal infra-red 

Interpretation of results obtainod in the 8-14 ym wavelengths 
has helped in outlining several interesting phenomena (Fig. 3.4, 
Going landward from coastal vjaters to the Dol marsh : 



’ Samples : 

1 

2D -2 mm 
in % 

2-0,2 mm 
in % 

200-20 m 
in % 

20-2 m 
in % 

inf. 2 m 
in % 

Organic 
matter 
in % 

Carbonate 
in % 

Water : 

content : 
in h : 

: 1 

- 

0,2 

27,8 

29,6 

42,6 

7,0 

27,9 

36,5 : 

: 2 

• - ” 

0,6 

42,5 

25,4 

31,5 

6,0 

30,2 

31,0 : 

: 3 

8,4 

31,4 

45,6 

5,7 

8,9 

5,5 

59,8 

18,5 ; 

: 4 

9,5 

39,0 

51,5 

- 

- 

1,5 

64,8 

12,0 : 

! 5 

10,7 

05,6 

3,7 

- 

. - 

0,5 

63,9 

1,0 : 

: 6 

- 

1.2 

60,5 

23,2 

15,1 

5,5 

42,5 

26,0 : 




Table • 3.4.E - Granulometry, organic matter content, carbonate and water contents for 
the various ground truth stations (of. Fig. 3.4. D] 
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Samples 

(A3 

(B3 

Stations 

(G3 

(D3 

(E3 

■(F3 

CoscinodiscuB lineatus 

+ 

+ 

+ 

■ + 



Navicula cancellata 

+ 



+ 

+ 


Gyrosigma fasciola 

•I* 

+ 

+ 

- 

+ 


Amphlprora sp. 


+ 

+ 

+ 

+ 

- 

CoBcinodiscus nitidus 

*!• 

+ 

4* 

~ 


- 

Pleurosigma angulatum 


■h 

+ 


•f- 

- 

Blddulphla regia 

- 

- 


** 

- 

- 

Mslosira dubia 




+ 

4* 



Table 3i4iF - Distribution of most abundant (+] diatomee species at 
various stations Cl] (cf. Fig 3. 4. 03 


(1) Measurements were made by Mr RICARO at the Laboratoire de 
Cryptogamie of Museum National d’Histolre Naturelle. 




Fig. 3.4.G. - Thermal infra-red mapping of the inter-tidal 
- zone off the Dol marsh CFig. 3.4.D.3 
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- arrival at the beginning of flavv tide of cold and clear 
waters in turbulent mixing with warmer bay waters which are 
heavily loaded with sediments. This contact takes place 
around the Chatelier and Rimains islands. 

~ digitations of warm water coming from the tidal flat are 
abundant at the limit between infra and inter-tidal zones. 

This flow takes its origin in the fish trap which is near 
the east bank of the Canal des Allemands channel. Orientation 
of this flow is easily explained by the fact that tidal 
currents were running N.V/ during the experience. 

- a temperature difference between warm diatomee covered 
"tangue" and colder barren flat. 

- a specific thermal signature for the mussel shells nurseries 
Cbouchots) which are higher in temperature than the diatomee 
covered "tangue" on w'hich they are settled 

- a thermal assymmetry between each side of the "V" shaped 

fish traps s the westward side being warmer than the east ward 

^ finally sand levees and barren fields of the Dol marsh are 
similar in surface temperature which- is much higher than 
for the other targets. 

b - 2 Visible and near infra-red 


The analysis of spectral signatures present in a test area of 
350 scan lines of 180 sample points has led to the study of a 
complete disjunctive coded data set Cc-fi- 2.4] with I = 63000 
sample points and J = 64 spectral coded equivalents. This matrix 
was analysed by the FRACAH automatic classification system. 

Interpretation of a classification in nine categories lead us to 
the following conclusions : 

- category limits are correct but some misclassif ications 
subsist in detailed category assignments 

- "tangue” is well discriminated by classification and 
subdivided in several classes depending on its water content 

- sandy levees as well as locally sand rich "tangue" masses 
are correctly outlined 

- diatomee covered areas are separated from barren ones but 
some misclassif ication occurs with cultivated fields of the 
Dol marsh 

- finally salt marshes and other parts of cultivated fields 

fall into the some category. > 
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4 - Criticism 


In conclusion, two points must be emphasized ; 

- correction for systematic bias in the visible and near infra-red 
was estimated on the tidal flat zone. Unfortunately it appeared 
that forward scattering of solar light was function of target 
albedo and surface texture so that trend correction was not 
complete for sea water and continental zones. It is this phenomenon 
which is responsible for most of the misclassif ications 

- however results were sufficiently good to assess the feasibility 
of mapping tidal coasts by automatic classification of remotely 
sensed data. In the present study thermal infra-red was found 

to be best for Knowledge of water movement and mixing in the west 
part of the Mont Saint-Michel bay while visible and near infra-red 
was better used for mapping tidal flats and associated organic 
films. 
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3. 5 - INVESTIGATION OF A COASTAL UPWELLING USING MOAA DATA 
5.1 - Data description 


Since 1972 NASA has regularly launched on a yearly basis a 
series of meteorological satellites which are run under the 
supervision of the National Oceanic and Atmospheric Administrati 
The NOAA satellites (currently N0AA5] are equipped with various 
electro optical systems. Among them the two Very High Resolution 
Radiometers (VHRR) are of interest for oceanographic application 

Satellite orbit is sun synchronous with an inclination of 101° 
over the equator ; this allows any point on earth to be surveyed 
twice a day. Time of overpass is approximately 9 AM and 9PM 
Universal Time for the coast of France. The geometry of this 
quasi-polar orbit tends to simplify the problem of geometric 
correction. 

The two radiometers operate in different spectral bands : 

- in the infra-red (IR) domain, from 10. 5y to 12y for the VH 

- in the visible domain, from OiSyto 0.7y for the VHRR2 thei 
identical systems consists of a 12.7 cm aperture elliptical 
mirror with a surface of 100.4 cm2. Rotation speed of the 
mirror is stabilized at 1 revolution per 150 milliseconds. 

The relative flux is reflected by the mirror and focused on a 
detector thus leading to a 0.6 str instantaneous field of view 
(IFOV). With the satellite cruising at an altitude of about 
'1460 k.ms, the ground resolution is n.'Sa km. Rotation speed of 
the mirror and trajectory velocity are synchronised so that no 
overlapping occurs in the relative geometry of consecutive scan 
lines. 

The horizon to horizon scanning is done through a 70° rotation 
angle, and during the remaining part of the 360° complete 
revolution, the mirrors scan various internal calibration 
equipement such as two black bodies maintained at a fixed 
temperature and a point in space opposite the sun. About twenty 
different calibration signals are acquired di/ring one revolution 
this information is used to convert the voltage output of the 
system into black body radiance. 

■ The signals transmitted by the satellite are acquired at the 
Meteorologie Nationals station (GEMS) in Lannlon (Bretagne). 

It is then digitized and preprocessed for data cleaning and 
compression of the calibration information. The data is avallab) 
in 8 bit coded digital form on magnetic tapes which are sent 
upon request to various users outside the met ecro logical communit 

At Ecole des Mines the main point of interest is deriving sea- 
surface temperatures maps at scales ranging from 1/100 000 to 
1/250 000 for various areas along the Atlantic and Mediteranean 
coasts of France. 
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A data file corresponding to an acquired image of the Earth 
Includes one record for each scan line. Each record in physically 
built as follows ; 


- 20 calibration ualues CIS bits] 

- 2068 pixel values [8 bits] for the IR channel 

- 2068 pixel values (8 bits] for the visible channel 

Processing is currently done on 800 to 2400 scan lines files. 
The number of lines for a complete coverage of France is 
approximately 1200. 

3. 5.2 - Data processing 


In the production of surface temperature maps, the first step 
is to convert the 8 bit digital counts into black body radiances. 
This is done using a locally linear approximation of the Boltzann 
formula based on the calibration information. From this step 
onwards, the data of interest, usually a subset of the ’entire 
scene, is handled in degree Clesius stored in binary -floating 
points form. 

The maps are produced as false-color displays or contour plots. 
The direct interpretation of these results has been very difficul 
due to the presence of spurious disturbances in the data. These 
problems were not experienced In other oceanographic applications 
such as the mapping and tracking of Gulf Stream meanders done by 
various NOAA experiment teams. 

This is due to the fact that in Europe' our scientific interests 
are more localized geographycally' and that we are looking at 
temperature gradients of a few degrees using the VHRR scanner 
which has a nominal performance between 1® C and 0.5° C. 

Due to these differences in scientific approach,- the oceanographi 
community in France Is much more demanding on the refinement of 
the data processing necessary to bring out useful details in the 
sea-surface temperature structure. 

Two types of disturbances have been identifleld 

- a background noise caused by transmission or analog to 
digital conversion 

- a striping eFfect due to sensor gain losses affecting the 
calibration signals. 

The striping effect Is similar in nature to the one currently 
experienced in LANDSAT imagery. However, in this case, one does 
not observe the six hand repetitlvity which is particular to the 
LANDSAT scanner technology. The jitter in average value from one 
line to the next displays on the contrary a random structure with 
characteristics varying in time during data acquisition. The 
variability can be as wide as it can be a purely random noise or 
a rather smooth modulation expanding over 10 to 20 scan lines. 

The lack of a typical frequency in this jitter led us to discard 
■frequency domain filtering as used in various oceanographic 
applications of LANDSAT imagery. 



other methods were used such as smoothing averages or modes of 
the scan line histogram time-series. These approaches did not 
prove sufficient to eliminate the scan line noise in NOAA data. 
This is why a new method was developed vjhich takes into account 
the complete structure of the marginal distribution of each soar 
line instead of only using a condensed information extracted fro 
this distribution such as an average or a mode veluev 

The so called background noise is essentially isotropic in natur 
and evenly- distributed in the image plane. It is responsible for 
the "salt and pepper" aspect of the display. Its elimination 
involves a straightforward two-dimensional weighted smoothing ■ 
algorithm. 

The cleaning effect is satisfactory especially when the algorith 
Is iterated several times, leading to a contouring effect in the 
display which is rather useful for oceanographic applications. 
The drawback of this approach is that the contrast at the bounds 
between the earth and the sea is smeared at each iteration with 
a consequent degradation of the geographic location of the 
coastline. 

This annoying effect Is of particular importance for coastal 
applications where ons would like to retain the original spatial 
definition of the coastline’ as produced by the ground resolutior 
of the NOAA scanner data C approximately 1 Kml . In order to meet 
this goal the smoothing was conditionned to the local gradient i 
relation with average RMS of the background noise. 

The unstriping and the iterative conditionned smoothing are 
sequentially applied to the data set of sea-sur-pace temperature; 
the final product is a smooth isocontour temperature field whicl 
can be easily interpreted by the oceanographers and have been 
found of great value when stored and compared for different . 
satellite overpass dates. 

5.3 - Significant results 

A - Interpretaion of imagery 

The data we have used for our interpretation is N0AA4 image? 
acquired on 30 July 1975 at 9h 17 mn TU. The generated maps 
(FigV 3. 5. A and B] cover the oceanic waters .surrounding the 
Finistere shoreline. In the visible band the gray scales ar? 
related to various radiance levels expressed in FT.L. In th; 
thermal infra-red band the gray scales are associated to 
radiometric temperatures measured in Celsius degrees (°C). 

The land features are characterized by radiances greater 
. than 10 FT.'L and temperatures greater than 2,'T°C, they appea: 
respectively in white [VIS] and black [IR] on the imagery^ 

In the visible band, some confusion is possible between lant 
and clouds and the coastline is rather fuzzy around the 
Crozon Cape. 
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Fig. 3.5.1 - Display of the NDAA-VHRR visible band 
Finistere coastal area. 
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Fig. 3.5.B - Display of NOAA-VHRR infra-red band 
Finistere coastal area. 
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The clouds with high radiance are easily outlined in the 
visible band Month of Finistere, organized in North-West 
strips. The corresponding features in the infra-red have 
temperatures under 13®C. However they are generally broader 
masses than in the visible, this is due to the inertia of the 
infra-red sensor whis is less capable of following instant 
contrasts. The resulting contamination has lead us to 
eliminate from our oceanographic interpreration most of the 
English channel. North of Finistere. 

Coastal fog is frequent along these coasts in summer, and 
their detection is more complex than that of the preceeding 
features. For example, in the Saint-Brieuc bay they could 
be interpreted in the visible as patches of turbid water 
whlth rather high radiances. This confusion is easily avoided 
when using the I.R because their apparent temperature is much 
lower than the average seasonal temperature of the .coastal 
waters of Brittany (typically 16°C). Probable presence of a 
small cloud of fog can be suspected North-VJsst of Ouessant. 

A cold surface front expands in 'the ocean waters- West of 
Finistere. Organized in a crescent like shapV,_frbm. the Raz 
point to the North cf Ouessant it is characterized surface 
temperatures between 13®C and 16®C. They are easily separated 
from the waters of the mouth of the English channel 
(> IG^C) and also from the coastal waters of South Brittany 
which are even warmer (above 17.5°C). Warm water pools are 
also locally trapped along the Finistere coast. South of the 
Raz point, in the Douarnenez bay and in the Brest sound. 

b - Ground truth checking 


Simultaneously with satellite overpass, a team of English 
oceanographers (1j was performing a detailed study of this 
cold front using classical surveying methods. The situation of 
the thermal gradients drawn from their maps is in excellent 
agreement with the N0AA4 observations. As far. as the numerical 
values are concerned, the difference between in-situ measuremenrd 
and satellite observations are no greater than 1°C. The radio- 
metric temperatures are always lower than sea-truth which could 
be due to atmospheric losses. 

Absolute temperature errors in the mapping of oceans thermal 
gradients using VHRR data are mainly due to the water content 
of the lower levels of the atmosphere and they can be corrected 
by the use of atmospheric models based for example on radio- 
sounding results. Large clouds are seldom an interpretation 
problem due to their straightforward separation in the visible 


(1) R.D. Pingree, P.R. Pugh, P.fl. Holligan and G.R. Forster 

Summer phytoplancton blooms and red tides along tidal fronts in the 
approaches to the English Channel. 

Nature - Vol 25B - PP B72-B77 



and In the infra-red but this is not the case for coastal 
fogs and scatters of small low clouds. The detection of 
these features is done by low clouds. The detection of these 
features is done by careful comparison between measured 
radiometric temperatures and seasonal average surface 
temperatures available for the surveyed area through 
classical physical oceanography observations. 

In the strong thermal gradient area of the coast of trinistere. 
measurements had also detected a fair increase in 'chlorophyll 
A concentration due to phytoplancton activity. These concen- 
trations should have modified the colour of the ocean by it 
is disappointing that the visible VHRR display does not show 
such a feacture. 

Several reasons can be conjectured : 

- the visible VHRR band, identical to the LANDSAT MSS5, 
is not suited to chlorophyll mapping in water ; 

“ light penetration (os 2 m) could blurr the precise 
definition of the colour gradients ; 

- the resolution (1 km] is too coarse regarding the 
detection of phytoplancton patches which follow the 
edge of the thermal front and are therefore of a v/idth 
comparable in size with the VHRR resolution. 

The results presented here show the utility of infra-red 
thermal imagery for oceanographic applications but they also 
suggest caution in the interp'retation of rather coarse 
resolution future satellites such as NIMBUS-G. 


3. 5.4 Conclusion 


Although the processing techniques outlined In this work, appear 
satisfactory for the available data some remarks can be made 
regarding the suitability of the resulting output for the user 
community's needs. 

The accuracy of sea surface temperatures obtained from NOAA data, 
although useful, is not up to the standards currently applied by 
oceanographers who are able to measure sea temperature to '1/100'th 
of a degree Celsius, using thermistors. The pf g degrees 

Celsius forecasted for future satellites will be closer to 
oQeanpgraphyc needs. 

Nonetheless thermal infra-red applications in general are more 
adaptated to 'oceanographic applications than, for Instance, Lendse 
data which provide ocean colour measurement at the surface or at e 
depth which is not always easily identified. 

On the other hand, the spatial (1 km square) resolution of the NOA 
VHRR is adequate for oceanographic applications, whilst Landsat de 
has the minimal resolution for coastal geomorphology applications. 



133 - 


The density of observations is almost excessive in oceanography 
and provides a great deal o-f redundant information. 

Finally, repetivity of coverage (twice daily for NOAA] provides 
time series wlch can be used for pratical monitoring of fast- 
changing situations. At the present time we have mostly processed 
day-time coverage (9 a.m.) but we are planning to use both day-time 
and' evening (3 p.m.3 data in conjunction with IHCMM data (12 a.m, - 
12 p.m.) for an evaluation of high frequency coverage which should 
enable us to monitor roost of the tidal cycle over a short period. 

In this last case, the low frequency or Landsat coverage makes it 
practically impossible to obtain a complete series since some tidal 
situations have not occured concurrenty with a Landsat overpass 
since we began to receive Landsat data. 
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